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ABSTRA CT 

A  digital  computer  program  is  described  which  simulates  the 
radar-target  engagement  providing  a  representation  of  the  detection, 
acquisition,  and  tracking  processes.  The  program  is  arranged  as  a 
time  simulation  of  the  engagement  between  a  radar  and  target,  taking 
into  account  the  detailed  characteristics  of  the  target  cross  section, 
radar  and  target  motion  throughout  the  engagement,  surface  clutter, 
atmospheric  attenuation,  and  radar  losses.  In  the  output  the  program 
provides  the  user  with  target  detection  probabilities  in  the  presence  of 
surface  clutter  as  well  as  receiver  noise,  radar  search  and  track  ac¬ 
curacies,  signal- to -noise  ratios,  target  characteristics  versus  time, 
angular  and  range  rates,  etc.  The  input  requirements  to  the  program 
are:  (1)  a  deck  of  parameter  cards  describing  the  radar  parameters, 
the  clutter  environment,  and  the  initial  radar- target  geometry,  (2)  a 
deck  of  cards  describing  the  target  motion  throughout  the  engagement, 
and  (3)  a  deck  of  cards  describing  the  target's  cross  section  versus 
aspect  angle.  Many  simplifications  to  the  inputs  are  allowed  for 
studying  and  isolating  various  parts  of  the  radar  problem. 
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PREFACE 

T.  is  paper  describes  a  digital  computer  program  which 
is  intended  to  serve  as  an  aid  in  the  design,  analysis,  and  evaluation 
of  radar  systems.  Techniques  for  analyzing  the  performance  of  a 
radar  without  the  aid  01  a  digital  computer  are  well  established 
and  provide  sufficiently  accurate  results  if  the  problem  is  not 
too  involved.  However,  the  problems  are  often  quite  complicated 
which  forces  the  radar  analyst  to  reduce  the  complexity  with 
simplifying  assumptions,  such  as,  specifying  the  target's  average 
cross  section  as  a  constant,  ignoring  the  effects  of  ati,  ^spheric 
attenuation  and  ground  clutter,  estimating  the  detection  probability, 
etc.  If  numerous,  accurate,  and  detailed  analyses  are  required  a 
digital  computer  program  must  be  used. 

The  purpose  of  this  program  is  to  simulate  the  radar- 
target  engagement  in  the  real  world  in  order  to  provide  a  representa¬ 
tion  of  the  detection,  acquisition,  and  tracking  processes.  In  the 
simulation  process  it  takes  into  account  the  detailed  characteristics 
of  the  target  cross  section,  radar  and  target  motion  through  the 
engagement,  surface  and  rain  clutter,  atmospheric  attenuation,  and 
radar  losses.  In  the  output  the  program  provides  the  user  with 
target  detection  probabilities  in  the  presence  of  surface  and/or 
rain  clutter  as  well  as  receiver  noise,  radar  search  and  track 
accuracies,  signal-to-noise  ratios,  target  characteristics  vs 
time,  angular  and  range  rates,  etc. 
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The  program  is  made  as  general  purpose  as  poss-ble  bv 
segmenting  it  into  subprograms  and  subroutines,  and  by  the  liberal 
use  of  data  inputs.  In  its  complete  form  broad  examples  of  its 
use  are  as  follows:  analyzing  the  performance  of  a  radar  against 
different  kinds  of  targets  for  different  sea  states;  optimizing 
the  yield  of  data  from  a  live  test  by  simulating  the  test  before¬ 
hand;  optimizing  the  scanning  pattern  or  any  other  parameter  of  a 
search  radar  in  various  situations;  and  providing  assistance  in 
determining  optimum  hardware  parameters  during  the  development  of 
a  radar  system.  If  a  complete  simulation  is  not  desired  the  target 
inputs  can  be  simplified  to  allow  the  radar  characteristics  to  be 
studied  separately. 

The  program  is  written  in  FORTRAN  II  for  the  IBM  7094 
computer.  It  compiles  in  approximately  two-tenths  of  an  hour 
and  computes,  for  an  average  number  of  cases,  in  two -hundredths 

of  an  hour. 

I  am  indebted  to  V.  Schwab,  G.  T.  Trotter,  and 
E.  Shotland  whose  contributions  represent  a  significant  part  of 


the  computer  program. 
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I.  INTRODUCTION 


The  digital  computer  program  described  in  this  report 
was  developed  to  calculate  the  performance  of  a  radar  against  a 
single  target  in  a  dynamic  situation.  The  program  is  hereinafter 
called  the  Radar  Analysis  Program  for  purpose  of  identification. 

1. l  The  Radar  Analysis  Problem 

The  radar  problem  that  is  to  be  analyzed  is  illustrated 
graphically  in  Figure  lrl.  Certain  complicating  factors  such  as 
time  variation  of  the  target  and  *radar  positions,  clutter  echoes, 
antenna  sidelobes,  and  changing  situations  prompts  the  radar 
analyst  to  use  a  digital  computer. 

Time  Variation.  The  prime  complicating  characteristic 
of  the  radar  analysis  problem  is  the  variation  of  the  positions  of 
radar  and  target  with  time.  If  the  target  and  radar  are  allowed 
to  move  throughout  an  engagement  then  a  detection  probability 
calculation  would  have  to  be  performed  for  every  increment  of 
time.  This  would  be  necessary  to  account  for  changes  in  the 
target's  cross  section,  range  to  the  target,  size  of  the  clutter 
echoes,  and  antenna  position.  The  change  of  target  cross  section 
as  a  function  of  time  is  calculated  with  a  fair  amount  of  precision. 
The  target  cross  section  as  a  function  of  aspect  angle  is  fed  into 
the  program  as  an  input.  The  actual  aspect  angle  as  a  function  of 
time,  as  calculated  during  the  engagement,  is  used  to  find  the 
corresponding  target  cross  section  by  table -look-up. 
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Fig.  1-1  RADAR  ANALYSIS  PROBLEM  IN  DYNAMIC  SITUATION 
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Surface  Clutter.  The  second  most  complicating  factor  in 
radar  performance  calculations  is  clutter,,  the  reflections  from  the 
background  land  or  sea.  The  magnitude  of  the  received  clutter  power 
is  calculated  each  increment  of  time  and  is  added  to  the  receiver 
noise  power.  This  total  noise  power  is  used  to  form  the  signal-to- 
noise-plus-clutter-ratio  which  is  the  prime  quantity  for  calculating 
the  detection  probability. 

The  most  popular  detection  theory  used  for  the  detection 
probability  calculation  was  developed  by  Marcum  and  Swerling.1 
However  this  theory  assumes  that  the  noise  is  random  from  pulse  - 
to-pulse  like  receiver  noise.  If  the  total  noi^e  includes  clutter 
echoes  the  Marcum-Swerling  theory  would  yield  an  erroneous  prob¬ 
ability  of  detection  since  clutter  echoes  are  in  general  not  randomly 
distributed  from  pulse-to-pulse .  To  improve  the  calculation  of  the 
detection  probability  for  a  target  in  a  clutter  background  a  new 
theory  is  used  which  treats  the  clutter  echoes  realistically,  i.e., 

some  degree  of  correlation  from  pulse-to-pulse.3  Both  detection 

j 

theories  are  incorporated  as  part  of  the  Radar  Analysis  Program; 
their  use  will  be  described  later. 

Antenna  Side  lobes.  The  clutter  problem  is  complicated 
further  when  the  antenna  sidelobes  are  included  in  the  simulation. 

For  example  under  situations  when  the  main  beam  of  the  antenna  is 


J.I.  Marcum  and  P.  Swerling,  "Studies  of  Target  Detection  by 
Pulsed  Radar,"  IRE  Transactions  on  Information  Theory.  Vol,  IT-6 
(April ,  1960) . 

8E.  Shot land,  "False  Alarm  Probabilities  for  Receiver  Noise  and 
Sea  Clutter,"  JHU/APL  Internal  Memorandum  BBD-1387 ,  October,  1964. 

-  3  - 
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pointing  at  a  target  at  a  high  altitude  the  clutter  echoes  are  not 
received  via  the  main  beam  but  with  the  sidelobes.  In  the  present 
simulation  the  sidelobes  are  considered  constant  at  a  level  specified 
in  the  program  input. 

Other  facets  of  the  radar-target  engagement  shown  in 
Figure  1-1,  or  otherwise  incorporated  in  the  Radar  Analysis 
Program,  will  be  described  in  the  remainder  of  the  report. 

Performance  Calculations.  Even  though  the  detection 
process  was  emphasized  in  the  above  it  is  not  the  only  performance 
characteristic  that  should  be  considered.  Other  performance 
characteristics  that  are  included  in  the  program  are  search  and 
track  accuracies  (in  range,  doppler,  and  angle),  range  resolution, 
and  target  identification  times. 

All  calculations  are  made  as  a  function  of  time  of  the 
engagement  at  any  desired  time  interval.  Any  number  of  runs  can 
be  performed  in  one  program  deck  set-up  to  analyze  radar  performance , 
or  target  characteristics,  as  various  parameters  (total  of  43)  are 
varied  one  at  a  time  or  together.  Examples  of  parameters  that 
are  often  varied  are  clutter  reflectivity  (or  sea  state) ,  radiated 
power,  antenna  gain,  sidelobe  level,  time  between  false  alarms, 
target  range,  and  angle  of  the  target's  plane  of  motion. 

1.2  Summary 

The  remaining  chapters  in  this  report  are  arranged 
according  to  the  major  divisions  of  the  computer  program,  namely, 
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inputs,  processing,  and  outputs.  The  Radar  Analysis  Program  is 
shown  in  block  diagram  form  in  Figure  1-2.  The  input  parameters 
are  arranged  in  five  groups  as  shown  and  are  discussed  in 
Chapter  II.  These  inputs  are  fed  into  the  processing  part  of  the 
program,  which  is  discussed  in  Chapter  III,  resulting  in  the 
program  outputs,  which  are  described  in  Chapter  IV. 

In  Chapter  II  the  input  parameters  to  the  Radar  Analysis 
Program  are  listed,  described,  and  supported  by  examples.  The 
prime  problem  that  will  be  used  as  an  example  throughout  Chapter  II 
and  the  remaining  chapters  is  the  detection  and  track  of  a  surface- 
launched  ballistic  missile  by  a  radar  mounted  in  an  aircraft.  This 
problem  has  all  the  elements  of  complexity  such  as  rapid  target 
motion,  fluctuating  target  cross  section,  high  clutter  background, 
etc.  and  is  the  best  example  for  illustrating  what  the  program  can 
and  cannot  do.  The  changes  required  on  the  input  cards  to  analyze 
other  problems  will  also  be  indicated.  The  inputs  required  for  the 
target  motion  and  the  target  cross  section  are  the  most  involved  and 
are  arranged  in  decks  (or  tables) ;  these  are  "looked-up"  in  the 
processing  operation.  The  other  parameters,  which  describe 
initial  radar-to-target  geometry,  the  clutter  and  noise  parameters, 
and  the  radar  characteristics,  are  fed  in  on  one  input  card  to  a 
parameter  for  a  total  of  43. 

In  Chapter  III  a  description  of  the  processing  and 
simulation  techniques  are  given.  An  overall  flow  diagram  of  the 
program  is  given,  followed  by  description  of  the  more  important 
subroutines  such  as  the  target  and  radar  motion  simulation,  the 
radar  simulation,  and  the  detection  routines. 

-  5  - 
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INPUT  (Chapter  II) 

r - 1 


Figure  1-2 

BLOCK  DIAGRAM  OF  THE  RADAR  ANALYSIS  PROGRAM 


TH*  OHMS  MORKIND  UNIVIWITY 

APPLIEC  PHYSICS  LABORATORY 

flLVt*  DANIN*.  MARYLAND 


In  Chapter  IV  the  output  parameters  (66  are  printed  out) 
are  described  followed  by  the  printed  and  plotted  output  for  the 
examples  given  in  Chapter  II.  Examples  will  be  shown  for  the  full- 
scale  dynamic  situation  and  the  simplified  simulation  in  which  the 
target  motion  is  static  and  the  target  r  section  is  constant. 

In  Chapter  V  the  significant  developments  of  this  study 
will  be  reviewed  along  with  the  limitations  of  the  existing  program 

and  areas  for  future  development. 

The  complete  program,  written  in  FORTRAN  II,  is  given 
in  Appendix  B.  A  literature  review  of  other  digital  computer 
programs  for  radar  analysis  or  simulation  is  given  in  Appendix  A. 


-  7 
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II.  PROGRAM  INPUTS 

The  program  data  inputs  consist  of  43  input -parameter 
cards  and  two  decks.  The  individual  parameter  cards  describe  the 
initial  radar-to -target  geometry,  the  radar  characteristics,  the 
clutter  and  noise  parameters,  and  program  control.  The  two  decks 
of  cards  form  two  tables  describing  the  target  motion  and  the 
target  radar-cross-section.  These  inputs  will  be  defined  and 
illustrated  by  examples.  The  prime  example  used  will  be  for  a 
radar  mounted  in  an  aircraft  with  a  surface-launched  missile  as 
the  target.  The  changes  required  in  the  inputs  for  other  examples 
will  be  indicated. 

•  The  last  section  in  this  chapter  will  describe  the 
input  parameters  and  procedures  used  to  control  the  type  of  output 
desired.  Two  types  are  available:  the  time  simulation  mode,  in 
which  the  radar  performance  is  analyzed  against  a  moving  target, 
and  the  detection-probability  range -prof ile  mode,  in  which  the 
target  motion  and  cross  section  are  unchanging. 

2. 1  Target  Motion 

The  target  motion  within  a  vertical  plane  is  described 
with  a  deck  of  input  cards.  On  each  card  in  the  deck  there  are 
eight  numbers  as  listed  and  defined  in  Table  2-1.  The  time  value 
referred  to  in  Table  2-1  may  have  a  zero  or  negative  value  on  the 
first  card  and  increase  from  card  to  card  with  any  desired  increment 
between  1  and  9999  seconds.  The  card  preceding  the  target  motion 
deck  has  the  time  increment,  total  number  of  cards,  and  the  name 

-  8  - 
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Card  number 

Time  (T) 

Down- range 
position 

Down  range 
speed 

Down  range 
acceleration 


Altitude  position 


Altitude  speed 


Table  2-1 

TARGET  MOTION  TABLE 


Maximum  of  322 


Time  (from  zero)  in  seconds 


Horizontal  distance  s.n  feet  of 
target  in  target  plane  at  time 
(T)  measured  from  initial 
target  position. 


Horizontal  speed  of  target  at 
time  (T)  in  feet  per  second 


Horizontal  acceleration  of  target 
at  time  (T)  in  feet  per  second 
per  second 


Target  altitude  (or  vertical 
position  in  target  plane) 
at  time  (T)  in  feet 


Vertical  velocity  of  target  at 
time  (T)  in  feet  per  second 


Altitude  acceleration 


Vertical  acceleration  of  target 
at  time  (T)  in  feet  per  second 
per  second 
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of  the  target.  The  parameters  that  describe  the  initial  radar- 
target  geometry,  discussed  in  a  subsequent  section,  assume  their 
given  values  on  the  time  specified  on  the  first  card  in  the  target 
motion  table. 

The  target  motion  for  the  airborne-radar/missile-target 
problem  is  described  by  a  deck  of  41  cards  using  the  format  in 
Table  2-1.  The  actual  table  is  printed  in  the  first  part  of  the 
program  output  discussed  in  Chapter  IV  and  will  therefore  not  be 
repeated  here.  It  describes,  however,  40  seconds  of  the  flight 
of  a  ballistic  missile  that  reaches  an  altitude,  in  this  time,  of 
24,000  feet  and  a  horizontal  range,  from  the  launch  point,  of 
10,000  feet. 

The  values  for  targets  other  than  a  surface-launched 
ballistic  missile  are  easily  prepared.  For  a  target  that  does 
not  move  during  the  engagement,  the  simplest  case,  only  two  cards 
are  required:  the  first  at  time  zero  and  the  second  at  the  maximum 
duration  of  the  problem.  Both  cards  would  have  the  desired  range 
and  altitude  values  with  the  velocity  and  acceleration  values  equal 
to  zero.  On  the  other  extreme  a  rapidly  moving  target  could  be 
simulated  which  could  not  possibly  be  achieved  by  a  practical 
target. 

The  only  restriction  which  this  procedure  has  is  that 
the  target  motion  is  restricted  to  a  vertical  plane.  For  missile 
targets  whose  trajectories  are  often  in  a  plane  this  is  not  a  problem 
but  it  may  be  for  aircraft  targets.  For  example,  in  the  present 
version  of  the  program  an  aircraft  target  cannot  turn  or  bank.  It 
may,  however,  dive  or  rise  in  a  flight  plan  confined  to  a  vertical 
Plai*e  •  s.‘  -  10  - 
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2.2  Target  Radar  Cross  Section 

The  radar  cross  section  of  a  target  is  a  terra  that 

relates  the  power  density  reflected  from  a  target  to  the  incident 

power  density  from  the  radar  and  is  defined  as  follows.1 

Power  delivered  per  unit 
solid  angle  in  the  direction 

.  of  the  radar _ 

0  w  Power  per  unit  area  incident 
on  the  target 

The  value  of  the  radar  cross  section  varies  as  a  function  of  the 
orientation  of  the  target  with  respect  to  the  radar,  the  polari¬ 
zation  of  the  radar,  the  radar  wavelength,  and  the  conductivity  of 
the  target's  surface.  In  this  program  the  radar  cross  section  of 
the  target  at  the  desired  radar  wavelength  is  required  and  is  fed 
into  the  computer  on  a  deck  of  cards. 

General  Description.  The  deck  of  cards  contains  the 
radar  cross  section  of  the  target  as  a  function  of  the  aspect 
angle,  for  both  vertical  and  horizontal  polarization.  The  aspect 
angle  is  defined  as  the  angle  between  the  imaginary  line  connecting 
the  radar  to  the  target  and  the  center  line  through  the  target.  For 
both  missiles  and  aircraft  the  center  line  runs  from  nose  to  tail 
with  the  r?ose-on  aspect  angle  taken  as  zero.  The  range  of  aspect 
angles  covered  are  from  0  to  180  degrees.  For  symmetrical  targets 
the  cross  section  values  for  the  aspect  angles  between  180  and  360 
degrees  are  the  same  as  the  values  for  the  corresponding  angles 
between  0  to  180  degrees. 

XR.  S.  Berkowitz,  ed. ,  Modern  Radar  (New  York:  John  Wiley  and 
Sons,  1966),  pg.  549. 
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The  radar  cross  section  table  is  stored  in  the  computer 
and  is  used  by  the  processing  part  of  the  Radar  Analysis  Program 
every  time  a  cross  section  value  is  desired  for  a  particular 
aspect  angle.  This  will  be  discussed  farther  in  the  next  chapter. 

Airborne-radar/missile  -target  example .  The  radar  cross 
section  deck  for  the  airborne-radar/missile-target  example  is 
partially  shown  in  Table  2-2  indicating  the  format  arrangement. 

The  complete  table  is  given  in  Chapter  IV;  it  was  prepared 
especially  for  this  report  and  does  not  represent  the  cross  section 
characteristics  of  a  known  missile.  The  values  of  the  cross 
section  table  were  plotted  in  Figure  2-1  to  show  the  general 
characteristics  of  the  target.  The  fine  detail  of  the  cross 
section  is  not  shown  in  this  figure  because  it  was  plotted  in 
increments  of  three  degrees.  As  shown  in  Table  2-2  the  actual 
cross  section  data  are  fed  into  the  program  in  .] -degree  increments 
over  aspect  angles  from  0  to  180°  for' Vertical  and: horizontal 
polarizations. 

This  type  of  data  is  available  from  the  Radar  Target 
Scattering  Range  (RAT  SCAT)3,  Holloman  Air  Force  Base,  New  Mexico. 
At  RAT  SCAT  the  radar  cross  section  of  various  targets  and  target 
models  is  measured  for  any  frequency  from  150  to  12,000  MHz  at 
arbitrary  polarizations.  A  program  was  developed  to  edit  the  RAT 
SCAT  data  (on  punched  paper  tape  or  magnetic  tape)  and  convert  it 
to  a  deck  of  cards  for  direct  input  to  the  Radar  Analysis  Program. 


®H,  C.  Marlow,  et  al .  ."The  RAT  SCAT  Cross-Section  Facility," 
Proceedings  of  the  IEEE,  Vol.  53,  No.  8.  Special  Issue  on  Radar 
Reflectivity,  (August  1965),  pp.  946*954. 

-  12  - 


EXAMPLE  OF  RADAR  CROSS  SECTION  INPUT  CARDS 


The  Johm  Hopkins  University 

APPLIED  PHYSIC*  LASOSATOSY 

Silver  Spring,  Maryland 


-  13  - 


First  two  cards  in  radar  cross  section  deck 


THC  JOHN*  HOPKINS  JNIVtRSITr 

APPLIED  PHYSICS  LASORATORY 

siivtp  opninq  Maryland 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 
8ilv*r  Spring.  Maryland 


Other  Targets.  The  cross  section  of  targets  other  than 
missiles  are  described  in  the  same  way  with  certain  limitations. 

For  an  aircraft,  for  example,  the  cross  section  as  a  function  of 
azimuthal  aspect  angles  for  both  horizontal  and  vertical  polari¬ 
zation  is  not  enough  to  completely  describe  the  target.  This  is 
because  the  cross  section  is  different  for  different  angles  off 
the  plane  of  the  wings.  This  is  not  the  case  for  missiles 
(neglecting  the  effect  of  fins)  and  other  targets  that  are  basically 
symmetrical  about  their  longitudinal  axis. 

The  radar  cross  section  input  technique  for  simple 
targets  like  a  sphere,  whose  cross  section  is  constant  as  a 
function  of  aspect  angle,  is  accomplished  by  entering  the  same 
cross  section  value  for  zero  degrees  aspect  angle  and  for  180 
degrees.  The  interpolating  routine  associated  with  the  part  of 
the  program  that  "looks-up"  this  table  will  then  choose  this 
constant  cross  section  value  for  all  aspect  angles.  This  same 
technique  can  be  used  to  describe  a  target  whose  cross  section  is 
constant  (or  can  be  assumed  constant)  over  a  limited  azimuth 
sector  changing  from  sector  to  sector  as  desired. 

2.3  Input  Parameter  Cards 

Following  the  target  motion  and  the  target  cross  section 
decks  is  a  set  of  43  input  parameter  cards.  These  cards  contain 
the  values  of  the  parameters  that  describe  the  initial  radar-to- 
target  geometry,  the  characteristics  of  the  radar,  the  clutter  and 
noise  parameters,  and  the  program  options;  these  parameters  are 

4  f  /n  j  r  ,  t’  i  g.  •.  /-  i. 

described  in  the  following  sections. 
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The  parameters  are  arranged  one  to  a  card  which  includes 
the  parameter  number,  the  FORTRAN  symbol,  the  r/alue,  the  conversion 
factor,  and  the  definition.  All  cards  are  printed  in  the  first 
pa*f  of  the  output  as  shown  in  Chapter  IV.  Parameters  that  are 
given  a  value  of  zero  can  be  omitted  from  the  input  parameter 
deck\ 

2.3.1  Initial  Radar-to-Target  Geometry 

The  initial  geometry  input  parameters  specify  the 
distances  and  angles  between  the  radar  and  target  at  the  start 
of  the  engagement,  and  the  velocity  vector  for  the  radar  motion 
during  the  remainder  of  the  engagement.  The  geometry  is  shown  in 
Figure  2-2.  /The  parameters  shown  are  defined  in  Table  2-3  with 
values  for  the  airborne -radar/missile -target  example. 

As  indicated  in  Table  2—3  the  parameters  for  the  airborne— 
radar/missile -target  example  describe  a  radar  in  an  aircraft  that 
is  flying  at  an  altitude  of  30,000  feet  in  level  flight.  The 
target  motion  throughout  the  engagement  is  confined  to  a  plane 
which  is  positioned  at  an  azimuthal  angle  of  30°  with  the  reference. 

The  changes  required  in  the  geometric  values  of  the 
input  parameters  for  problems  other  than  the  airborne-radar/ 
missile -target  problem  are  apparent.  If  the  radar  were  on  a  ship, 
for  example,  the  altitude  ZLOIT  would  be  reduced  to  zero  and  the 
velocity  WA  reduced  to  say  70  feet  per  second.  The  other  angles 
and  distances  could  be  changed  as  desired.  Note,  however,  that 
in  all  cases  the  target  motion  (described  in  Section  2.1)  through- 

£  j.  '  .  c  ’( '  ° 

out  the  remainder  of  the  problem  is  confined  to  a  vertical  plane, 
which  is  not  a  limitation  for  most  problems. 
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Table  2-3 

INITIAL  RADAR-TO-TARGET  GEOMETRY 
INPUT  PARAMETERS 


PARAMETER 

DEFINITION 

EXAMPLE 

VALUE 

XLOIT 

Distance  between  radar  and  target 
along  azimuth  reference  line,  in 
nautical  miles 

30  nm 

YLOIT 

Distance  between  radar  and  target, 
perpendicular  to  azimuth  reference 
line,  in  nautical  miles 

30  nm 

ZLOIT 

Altitude  of  radar,  in  feet 

30,000  ft 

AZTAR 

Angle  between  target  plane  of 
motion  and  zero-azimuth  reference 
line,  in  degrees 

CO 

o 

0 

VVA 

Velocity  vector  of  radar  platform 
at  start  and  throughout  engage¬ 
ment  ,  in  feet  per  second 

700  fps 

ELVA 

Elevation  angle  of  VVA,  in  degrees 

0° 

AZVA 

Azimuth  angle  of  VVA,  in  degrees 

60° 

OLDK 

Switch,  0  or  1,  multiplying  VVA 

1 
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2.3.2  Radar  Characteristics 

The  radar  characteristics  are  fed  into  the  program  on 
input  cards,  one  parameter  to  a  card,  and  can  be  grouped  into 
search  radar  and  track  radar  parameters. 

Search  Radar.  The  search  radar  parameters  are  listed 
in  Table  2-4  with  descriptions  and  example  values.  These  values 
will  be  used  in  the  airborne-radar/missile-target  example  but  are 
otherwise  completely  arbitrary.  As  listed  these  parameters 
describe  a  5000  MHz  (C-band)  radar,  radiating  500  kilowatts  of 
power  in  pulses  that  are  10  microseconds  wide  and  occur  at  a 
repetition  rate  of  1000  per  second.  Pulse  compression  is  used 
with  a  time  bandwidth  product  of  100.  This  indicates  that  the 
10  microsecond  transmit  pulse  will  be  compressed  with  a  receiver 
that  is  matched  to  the  code  within  the  transmit  pulse,  to  a  pulse 
that  is  .1  microsecond  wide.  The  system  loss  factor  is  listed  as 
.1  (-10  dB)  and  includes  atmospheric  loss,  radar-line  loss,  beam- 
shape  loss,  and  scanning  loss.3 

The  remaining  search  parameters  describe  the  antenna 
and  its  scanning  characteristics.  A  vertically  polarized  antenna 
was  chosen  which  is  5.6-feet  square  (generating  a  2°  pencil  beam) 
with  an  average  sidelobe  level  of  30  dB  below  the  mainlobe  gain. 

The  antenna  beam  will  be  scanning  26,6°  sector  in  azimuth  and  a 
30°  sector  in  elevation,  in  one  secord.  The  time  to  scan  the 
sector,  HTSS,  is  the  program  clock  and  will  specify  the  increment 
of  time  at  which  the  radar  calculations  are  performed  by  the  program. 

aD.  K.  Barton,  Radar  System  Analysis  (Englewood  Cliffs:  Prentice- 
Hall,  Inc.,  1964) ,  p.  140. 
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Table  2-4 

SEARCH -RADAR  INPUT  PARAMETERS 


PARAMETER  DESCRIPTION  EXAMPLE 

BMC  Radar  frequency  in  megacycles  5000  MHz 

PPEAK  Peak  power  in  megawatts  .5  Mw 

TDWELL  Width  of  transmit  pulse  in  10  jis 

microseconds 

FR  Pulse  repetition  frequency  in  Hz  1000  Hz 

FNOISE  Receiver  noise  figure  3  (5  dB) 

SYSLF  System  loss  factor  .1 

PCN  Pulse  compression  ratio  100 

POL  Antenna  polarization;  .  one  for  1 

horizontal  and  0  for  vertical* 

HAPER  Horizontal  aperture  of  antenna  5.6  ft 

in  feet 

VAPER  Vertical  aperture  of  antenna  in  5.6  ft 

fdet 

AAPEFF  Antenna  aperture  efficiency  .65 

ETA  Si delobe  to  mainlobe  antenna  gain  .001 

ratio  in  decimal  form 

XIHDEL  Azimuth  search  sector  in  degrees  26.6° 

XIVDEL  Elevation  search  sector  in  degrees  30° 

HTSS  Time  to  scan  the  sector;  also  the  1  sec 


time  increment  between  all  calculations 


POL  -  0.  actually  indicates  that  the  antenna  polarization  vector 
is  parallel  with  the  plane  of  target  rotation,  and  POL  =  1  when 
the  polarization  vector  is  perpendicular  to  the  plane  of  rotation. 
Therefore  for  a  ballistic  missile  and  a  vertically  polarized  radar 
POL  -  0,  while  for  an  aircraft  target  and  a  vertically  polarized 
radar  POL  =  1. 
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The  search  radar  parameters  will  be  used  in  the  standard 
radar  equation  for  calculation  of  the  signal-to-noise  ratio.  The 
parameters  may  be  modified  to  calculate  the  performance  of  radars 
with  modulations  different  from  the  pulse  radar  depicted  here.? 

A 

For  example,  when  analyzing  a  CW  (continuous-wave)  radar  the  peak 
power  PPEAK  is  equal  to  the  average  power;  the  pulse  width  TDWELL 
is  equal  to  the  time  the  beam  remains  on  the  target;  and  the  pulse 
repetition  frequency  FR  is  adjusted  to  make  the  number  of  pulses 
that  hit  the  target  equal  to  one. 

Track  Radar.  The  track  radar  input  parameters  for 
the  Radar  Analysis  Program  are  listed  in  Table  2-5,  These 
parameters  are  used  to  calculate  the  angle  tracking  accuracies 
according  to  equations  by  Barton.®  The  parameters  are  discussed 
completely  in  Barton6  for  monopulse  and  conical  scan  trackers. 

The  values  of  the  parameters  indicated  in  Table  2-5  in  the  example 
column  are  for  a  monopulse  tracking  radar. 

2.3.3  Clutter  and  Noise  Parameters 

The  clutter  and  noise  input  parameters  are  required  to 
determine  the  method  for  calculating  the  amplitude  of  the  clutter 
and  the  detection  theory  to  be  used  in  calculating  the  probability 
of  detection.  The  six  input  parameters,  with  a  brief  description 
and  example  values  for  the  airborne -radar/missile -target  example, 

*J.  J.  Bussgang,  >  et  al., ,  "A  Unified  Analysis  of  Range  Performance 
of  CW,  Pulse,  and  Pulse  Doppler  Radar,"  Proceedings  of  the  IRE, 
Vol .  47,  (October  1959),  pp.  1753-1762. 

6 Barton,  o£.  cit. ,  p.  279. 

6  Ibid.  ,  pp.  263-315. 
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PARAMETER 

FNTRK 


CRL 

BN SR 

FKS 

B1 

ATETS 

ATEPS 

ATEOS 


Table  2-5 

TRACK-RADAR  INPUT  PARAMETERS 


DESCRIPTION  EXAMPLE 

Track  system  factor;  correcs  for  2 

detector  and  type  of  tracker; 

Monopulse  (coherent  det)  =  1 
Monopulse  (sq.  law  det)  =  2 
Conical  scan  (linear)  =  2 
Silent  lobing  (coherent)  =  2 
Silent  lobing  (sq.  law)  =  4 


Beam  cross-over  loss  factor;  equals  1. 

1  except  for  conical  scan  linear 
detector  in  which  case  it  equals  2. 

Nominal  noise  bandwidth  of  angular  1 

circuits;  /b 

n 

Normalized  slope  factor;  1,57  for  1.57 

monopulse 

Box  car  bandwidth,  cps;  PRF  for  1000. 

pulse  system;  bandwidth  of  doppler 
filter  for  pulse  doppler  system 

Angular  tracking  error  due  to  target  5 

motion,  in  milliradians 

Angular  tracking  error  due  to  .5 

platform  motion,  in  milliradians 


Angular  tracking  error  due  to  other 
causes,  in  milliradians 
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are  given  in  Table  2-6. 

The  quantity  aQ  listed  in  Table  2-6  is  defined  as  the 

normalized  clutter  cross  section  and  is  equal  to  the  radar  cross 

section  of  the  clutter  medium  per  unit  area  intercepted  by  the 

antenna  beam.  A  description  of  aQ  can  be  found  in  Sktilnik7 

which  includes  representative  values  for  different  surfaces,  A 

more  detailed  treatment  of  radar  backscatter  from  land,  sea,  or 

atmosphere  can  be  found  in  the  classified  literature.* 

The  selection  or  specification  of  the  clutter  and  noise 

parameters  for  any  given  case  is  best  described  with  the  aid  of 

Figure  2-3.  A  flow  chart  is  shown  with  a  box  for  each  of  the  six 

parameters;  the  double  lines  between  boxes  indicate  the  path  taken 

by  the  airborne -radar/missile -target  example.  The  time  between 

false  alarms  TAUFA  is  specified  first.  An  option  for  determining 

the  value  of  oq  is  then  reached.  Three  options  are  available: 

aQ  is  zero,  that  is,  there  is  no  clutter;  aQ  is  given,  in 

which  case  a  value  is  selected;  or  a  is  to  be  calculated 

o 

internally.  In  this  latter  case  the  sea  state  (if  in  fact  the 

operation  is  taking  place  over  the  sea)  is  specified  and  the  value 

of  a  is  calculated  with  an  empirical  equation  based  on  the  sea 
o 

state,  angle  of  incidence  of  the  radar  energy  to  the  surface  of 
the  sea,  and  the  radar  polarization. 

Chain  clutter,  v/hich  is  added  to  receiver  noise,  can  be 

'  .  V 

included  by  means  of  the  last  two  parameters  in  Table  2-6.) 

/ 

7M .  I.  Skolnik,  Introduction  to  Radar  Systems  (New  York:  McGraw- 
Hill  Book  Company,  Inc.,  1962),  p.  523, 

8  F.  E.  Nathanson,  ed. ,  "Report  of  Radar  Clutter  Signal  Processing 
Committee:  Part  I,  Radar  Clutter  Effects  (U) , "  TG  842-1,  The 
Johns  Hopkins  University  Applied  Physics  Laboratory,  September 
1966,  CONFIDENTIAL. 
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Table  2-6 


CLUTTER  AND  NOISE  INPUT  PARAMETERS 


PARAMETER 


DESCRIPTION 


EXAMPLE 


TAUFA 


False  alarm  time  in  seconds 


10  sec 


SIGOPT 


aQ  calculation  option: 

-1  »  no  clutter 
0  =  calculate  aQ 
1  =  given  as  SIGZ 


SIGZ 


FDMS 


CASE 


Value  of  a „  in  dB  when  SIGOPT  =>  1 
o 

Sea  state: 

1  =  rough  sea  (Beaufort  state  5, 
10  ft  wave  height) 

0  =  calm  sea  (Beaufort  .5, 

1  ft  wave  height) 

Frequency  diversity: 

1  «  Use  Marcum  and  Swerling 
0  «  Use  Shotland  theory 

Swerling  target  case  designation: 

0  «  nonfluctuating  target 

1  =»  scan-to-scan  fluctuation 

with  many  nulls 

2  =*  pulse -to-pulse  fluctuation 

with  many  nulls 

3  =  scan-to-scan  fluctuation 

with  a  few  nulls 

4  =  pulse -to -pulse  fluctuation 

with  few  nulls 


RAIN 


SUMSJG 


For  rain  clutter  at  the  target  RAIN  =  1. 
Otherwise  equals  zero. 

2 

Backscatter  coeff  icient  Zcr  in  meter  / 
meter^  in  decibels.  (See  Skolnik,  p.  539) 


wf 

VAV  ..  -4c 
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The  FDMS  option  is  now  reached  which  determines  whether 
the  Marcum-Swerling  theory  or  the  Shot  land  theory  is  to  be  used. 

In  essence  the  parameter  FDMS  is  not  an  option  bu*  a  radar  parameter 
that  describes  whether  frequency  diversity,  i.e.,  frequency  jumping 
from  pulse -to -pulse ,  is  used  or  not.  If  frequency  diversity  is 
used  then  the  clutter  is  considered  random  from  pulse -to-pulse9 
and  Marcum-Swerling  detection  theory  is  used.  If  the  radar  does 
not  use  frequency  diversity  then  the  clutter  is  not  random  from 
pulse -to -pui&e  which  means  that  the  Shot land  det  ction  theory  must 
be  used.  Note  that  if  there  is  no  clutter,  SIGOPT  =  -1,  then  the 
FDMS  option  is  bypassed  and  Marcum-Swerling  theory  is  used  in 
calculating  the  detection  probability. 

2.3.4  Program  Options 

There  are  two  options  for  operating  the  Radar  Analysis 
Program  which  are  selected  by  means  of  the  TARGOP  parameter  listed 
in  Table  2-7. 

Time  Simulation  (Primary  Option).  The  primary  option 
starts  from  TSTART  and  ends  at  TSTOP  and  runs  as  a  time  simulation 
of  the  radar-target  engagement.  Calculations  of  the  radar 
detection  probabilities  and  other  performance  characteristics, 
are  made  every  time  the  antenna  scans  past  the  target.  This 
mode  is  selected  when  the  option  TARGOP  is  equal  to  zero  or 
some  negative  value. 


* V.  W.  Pidgeon,  "Time,  Frequency,  and  Spatial  Correlation  of 
Radar  Sea  Return,"  A  Technical  Note  for  Use  of  Space  Systems 
for  Planetary  Geology  and  Geophysics,  The  American  Astronaut ical 
Society,  May  1967. 


i 

1 

] 

] 

1 

] 

] 

] 

] 

] 

] 


TMI  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 
Silvc*  Spring,  Maryland 


Table  2-7 

INPUT  PARAMETERS  FOR  PROGRAM  OPTIONS 


PARAMETER 


DESCRIPTION 


TSTART 

TSTOP 


TARGOP 


BXLOIT 

XLOIT 


Program  starting  time  in  seconds 
relative  to  zero  time  on  the 
target  motion  input 

Program  stopping  time  in  seconds 
relative  to  zero  time  on  the 
target  motion  input 

Program  option: 

-x  or  0  -  Time  simulation 
+R  -  Range  profiles 

where  R  is  the  range 
increment 

Maximum  range  value  for  performing 
calculations  in  the  range  profile 
option 

Minimum  range  value  for 
calculations  in  the  range  profile 
option 


EXAMPLE 
0  secs 


39  secs 


0 

(2  n.m.  in 
2nd  run) 


78  n.m. 

(in  2nd  run) 


30  n.m. 
(zero  n.m. 
in  2nd  run) 
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Ran^e  Profile  (Secondary  Option).  The  inputs  described 
to  this  point  are  those  required  for  the  time  simulation  option. 
However,  the  program  can  also  be  used  for  plotting  the  radar 
performance  versus  range  for  a  stationary,  constant —cross-section 
target  by  simply  specifying  TARGOP  as  +R.  For  a  stationary  target 
with  a  constant  cross  section  the  two  input  decks  are  greatly 
simplified:  the  target-motion  input  deck  is  reduced  to  one  card 

indicating  the  desired  altitude,  and  the  cross  section  input  deck 
is  reduced  to  one  card  indicating  the  desired  cross  section.  The 
program  then  calculates  the  detection  probability  and  other 
parameters  as  a  function  of  range  in  increments  of  R.  These 
calculations  versus  range  are  called  range  profiles  and  will  be 
discussed  further  in  Chapter  IV. 

Deck  Arrangement.  The  input -deck  arrangement  for 
running  the  example  problems  discussed  in  this  chapter  is  listed 
in  Table  2-8.  Two  runs  are  made:  (1)  the  airborne -radar/missile - 
target  example  in  a  full  time  simulation  mode,  and  (2)  the  airborne 
radar  versus  a  stationary  target  in  a  range  profile  run.  The 
items  listed  in  Table  2-8  were  discussed  in  this  chapter  except 
for  the  control  cards  and  cover  cards.  The  control  cards  include 
the  computer  job  cards  and  loading  cards.  The  cover  cards  are 
associated  with  the  data.  There  is  a  target-deck  cover  card 
which  indicates  the  target  name,  number  of  cards,  and  the  time 
increment.  The  cross  section  deck  has  both  cover  and  trailer 
cards.  The  former  indicating  the  name  of  the  target  and  the 
increment  of  the  aspect  angle;  the  latter  indicating  the  end  of 
tho  cross  section  deck. 
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Table  2-8 

ARRANGEMENT  OF  INPUT  DECKS 
FOR  MULTIPLE  RUNS 


Description  of  Approximate  Number 

Cards  of  Cards 

Control  cards  3 

Program 

Data  (Time  simulation  run) 

Target  Motion  Deck  42 

Cross  Section  Deck  300 

Input  Parameters  43 

Run  1  End  1 

Run  2  Leader  (Range  Profile  Run)  1 

Stationary  Target  Card  2 

Cross  Section  Card  3 

Input  Parameters  Different  from  Run  1  3 
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2.4  Summary 

The  input  cards  to  the  Radar  Analysis  Program  have  been 
described  in  some  detail.  The  target  motion  in  a  vertical  plane 
throughout  the  engagement  and  the  target's  radar  cross  section 
versus  aspect  angle  are  described  on  cards  to  practically  any 
degree  of  detail.  For  the  example  problem  studied  here  the 
target  motion  deck  comprised  of  41  cards  describes  a  missile 
launched  from  the  ground  and  reaching  an  altitude  of  24,000  feet 
in  40  seconds.  The  target's  radar  cross  section  is  described  on 
300  cards  with  values  every  .1  degree  for  aspect  angles  from  0 
to  180°  for  vertical  and  horizontal  polarization.  The  cross 
section  values  vary  over  50  dB  from  -25  dB  to  +25  dB  relative 
to  one  square  meter. 

The  43  input  parameter  cards  describe  the  initial 
positions  of  the  radar  and  target,  the  radar  parameters,  the 
clutter  and  noise  parameters,  and  the  program  mode  options.  For 
the  example  problem  these  parameters  describe  an  airplane  flying 
at  an  altitude  of  30,000  feet  away  from  a  missile  launched, 
essentially  straight  up,  at  an  initial  range  of  43  nautical  miles. 
The  radar  in  the  airplane  has  an  antenna  5.6  feet  on  a  side 
operating  at  a” frequency  of  5000  MHz,  radiating  500  kilowatts 
Of  power  in  a  10  jis  pulse.  The  radar  employs  pulse  compression 
and  frequency  jumping  from  pulse-to-pulse  to  reduce  the  effects 
of  clutter  contributed  by  a  rough  sea.  The  antenna  is  scanning 
a  sector  26.6°  in  azimuth  and  30°  in  elevation,  centered  on  the 
target,  in  one  second. 
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The  calculations  performed  on  these  inputs  will  be 
described  in  the  next  chapter.  The  resulting  outputs  are 
discussed  in  Chapter  IV.  / 
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III.  PROCESSING 

The  processing  part  of  the  Radar  Analysis  Program 
includes  those  calculations  and  analytical  models  required  to 
determine  target  cross  section  characteristics  versus  time, 
detection  probabilities,  time  required  for  detection,  tracking 
accuracies,  and  search  accuracies.  The  calculations  have  been 
divided  into  three  groups  for  purposes  of  discussion:  geometry 
calculations,  radar  calculations,  and  detection  calculations. 

These  three  blocks  of  calculations  are  shown  in  Figure  3-1,  the 
Radar-Analysis-Program  flow  diagram,  which  is  discussed  next. 

3.1  Flow  Diagram 

The  simplified  flow  diagram,  Figure  3-1,  indicates  the 
general  calculating  procedure  for  the  time  simulation  or  main  mode 
of  operation.  The  range-profile  mode  of  operation  is  a  special 
case  and  can  also  be  described  with  the  aid  of  Figure  3-1. 

Time  Simulation.  The  main  mode  of  operation  reads 
the  input  cards  described  in  Chapter  II,  stores  this  information, 
and  prints  it  out.  The  geometry  calculations  are  performed  next 
at  time  TSTART .  In  these  calculations  the  target-motion  input 
deck  and  the  initial  radar-to-target  geometry  parameters  are 
used  to  calculate  the  ranges,  angles,  and  range  rates  between 
the  radar,  target,  and  other  points  of  interest.  The  aspect 
angle  determined  in  this  set  of  calculations  allows  the  cross 
section  of  the  target  to  be  determined  by  looking  up  the  value 
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Start 


READ  INPUTS 

1. 

Target  motion  deck 

2. 

Target  radar-cross  section  deck 

3. 

Parameter  cards 

GEOMETRY  CALCULATIONS 

Aspect  angles 
Range 

Range  rates,  etc. 

Look  up  target  cross  section 


Read  new 
inputs 
and  repeat 


RADAR  CALCULATIONS 

Clutter  power 
Signal-to-noise  ratios 
Tracking  accuracies,  etc. 


Increase  time 
by  sector- 
scan  time 
and  repeat 


DETECTION  CALCULATIONS 

Marcum-Swerling 

or 

Shotland 


Is  this' 
;he  last, 
^^run?/ 


Is  time 
\TSTOP 

X?  / 


PRINT  OUTPUT 


Figure  3-1 

RADAR  ANALYSIS  PROGRAM  FLOW  DIAGRAM 
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in  the  radar-cross-section  input  table.  These  geometry  and  cross- 
section  values  are  used  in  the  next  set  of  calculations,  called 
radar  in  Figure  3-1,  to  determine  the  signal-to-noise  ratios  and 
the  tracking  accuracies  (to  .,ame  a  few) .  The  signal-to-noise 
ratio  and  the  clutter  and  noise  input  parameters  are  then  used 
in  the  detection  calculations  for  finding  the  single-scan  detection 
probabilities . 

The  time  step  of  this  set  of  calculations  is  then 
compared  with  the  last  time  for  calculations,  TSTOP.  If  TSTOP 
has  not  been  reached  the  time  variable  is  increased  by  the 
increment  of  time  required  to  scan  a  sector,  HTSS,  and  the 
calculations  are  repeated.  When  TSTOP  is  finally  reached  the 
output  is  printed  signifying  the  end  of  one  run.  Then  if  a  new 
set  of  input  parameters  exists  in  the  deck  they  are  read  and  the 
entire  processing  operation  is  repeated. 

Range  Profiles.  If  the  range  profile  option  is  selected, 
by  specifying  TARGOP  as  +R,  the  processing  operation  is  somewhat 
simplified  and  the  format  of  the  output  is  different.  Referring 
to  Figure  3-1,  the  input  is  read  as  before  except  now  the  target- 
motion  deck  and  the  cross-section  deck  consist  of  only  one  card 
each.  The  parameter  cards  are  essentially  the  same  as  for  the 
time  simulation  except  that  the  parameter  TARGOP  equals  +R. 

(Those  input  parameters  which  do  not  apply  to  the  range-profile 
case,  such  as  TSTART  and  TSTOP,  need  not  be  changed  or  omitted 
since  they  are  automatically  avoided  by  the  program.)  This  set 
of  inputs,  as  described  in  Chapter  II,  specify  a  target  whose 
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altitude  is  constant,  and  whose  cross  section  is  constant  and 
independent  of  aspect  angle. 

The  processing  for  the  range-profile  option  omits 
most  of  the  geometry  calculations  and  performs  the  radar  and 
detection  calculations  as  a  function  of  range.  Instead  of  the 
comparison  of  time  with  TSTOP  as  shown  in  Figure  3-1  the  range - 
profile  processing  compares  the  range  with  the  maximum  specified 
range  BXLOIT.  If  the  range  is  smaller  than  BXLOIT  it  is  incremented 
by  +R  and  the  calculations  are  repeated.  If  the  range  equals 
BXLOIT  the  run  is  printed  out  and  the  next  run  is  processed. 

Subroutines.  The  overall  processing  as  described  above 
is  controlled  by  the  main  program.  All  the  calculations  are 
performed  by  subroutines.  When  the  calculations  are  completed 
control  is  returned  to  the  main  program.  The  main  program  is 
then  used  for  printing  and  plotting  the  output  data. 

There  are  two  subroutines  that  are  associated  with  the 
main  program  that  do  not  perform  calculations.  These  are  TARGIN 
and  CROSIN  which  are  used  in  the  first  block  shown  in  Figure  3-1 
to  read  in  the  target-motion  and  cross-section  input  decks.  They 
also  are  responsible  for  printing  out  this  input  data  in  the  first 
part  of  the  output.  The  main  program  and  the  subroutines  are 
given  in  their  FORTRAN  II  language  in  the  Appendix  B.  The  first 
page  of  Appendix  B  lists  the  parts  of  the  program  in  the  order 
called. 
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The  actual  calculations  performed  will  be  discussed  in 
the  next  three  sections.  The  final  form  of  the  output  can  be 
found  in  Chapter  IV. 

3 . 2  Geometry  Calculations 

The  geometry  calculations  are  performed  by  the  subroutine 
GEOM  and  its  associated  subroutines:  TARGET,  AIRCFT ,  five  vector 
subroutines,  and  RATSCT .  (See  Appendix  B,  Section  2.0,  for  a 
listing  of  these  subroutines.) 

The  geometry  subroutine  is  called  by  the  main  program 
each  increment  of  time  to  solve  the  geometry  illustrated  in  Figure 
3-2.  The  geometry,  as  indicated,  is  solved  in  three  dimensions 

l 

over  a  spherical  earth. 

Four  points  in  space  are  involved  in  the  calculations: 
the  radar  location,  the  target  location,  the  point  of  reflection 
on  the  earth's  surface,  and  the  clutter  spot.  The  procedure 
Involved  in  calculating  the  pertinent  ranges,  range  rates,  and 
angles  between  these  points  is  as  follows.  The  target  position 
and  velocity  at  the  time  of  interest  are  gotten  from  the  target- 
motion  input  table,  using  interpolation  if  required;  this  is 
performed  by  subroutine  TARGET.  The  radar  position  at  the 
calculating  time  is  determined  by  the  subroutine  AIRCFT,  which 
determines  the  distance  traveled  due  to  the  input  velocity  WA 
specified  for  the  radar  platform.  Then  the  ranges,  range  rates, 
and  angles  connecting  the  radar  and  target  can  be  computed.  They 

*A  4/3  earth  radius  is  used  to  account  for  refraction. 
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are  computed  using  subroutines  for  each  of  the  vector  operations 
involved;  namely,  finding  the  cross  product  between  two  vectors 
(CROSS) ,  finding  the  dot  product  between  two  vectors  (DOT) ,  reducing 
a  vector  to  its  unit  vector  and  magnitude  (UNIT) ,  multiplying  a 
vector  by  a  scalar  (MULT) ,  finding  the  angle  of  a  vector  given 
the  sine  and  cosine  (ANGLE) ,  and  finding  the  elevation  and 
azimuth  angle  of  a  vector  with  reference  to  a  given  coordinate 
system  (TRIAD) . 

The  same  procedure  is  used  to  find  the  ranges,  range 
rates,  and  angles  between  the  radar  and  the  reflection  point,  and 
the  radar  and  the  clutter  spot.  The  reflection  point,  located 
at  the  surface  of  the  earth,  comes  into  play  when  the  radar's 
antenna  is  broad  enough  tc  permit  energy  to  be  reflected  off  the 
surface  of  the  earth  and  up  to  the  target.  Two  types  of  reflections 
are  considered:  single  oounce,  where  the  radar  energy  follows 
the  reflected  path  to  the  target  and  the  direct  path  back,  and 
double  bounce,  where  the  energy  travels  the  reflected  path  to 
and  from  the  target.  The  clutter  spot  is  simply  the  point 
on  the  surface  at  a  range  equal  to  the  range  to  the  target. 

The  aspect  angles,  also  shown  in  Figure  3-2,  are 
calculated  by  the  geometry  subroutine  and  are  used  directly  in 
finding  the  target's  cross  section  this  instant  of  time. 
Subroutine  RATSCT  is  used  to  look  up  the  cross  section  table. 

The  cross  section  value  selected  is  the  value  opposite  the  closest 
angle  to  the  desired  angle. 
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There  are  a  total  of  31  parameters  (many  with  three 
components)  which  are  calculated  and  stored  by  the  geometry 
subroutine.  However,  at  the  present  time  only  the  most  important 
ones  are  printed  as  described  in  Chapter  IV  on  output. 

3.3  Radar  Calculations 

The  radar  calculations  are  performed  each  time  step 
by  subroutine  DAVE  following  the  execution  of  the  geometry 
subroutine.  These  calculations  w'll  be  briefly  described  here 
with  the  details  being  available  in  Section  3.0  of  Appendix  B. 

Antenna.  The  physical  size  of  the  antenna,  the  efficiency 
and  the  radar  frequency  are  used  to  calculate  the  antenna  gain  and 
the  beamwidth.  The  value  for  the  beamwidth  is  then  used  1th  the 
given  values  for  the  search  sector  and  the  time  allowed  to  scan 
the  sector  to  determine  the  scan  rate.  The  scan  rate  is  used 
with  the  given  value  for  the  pulse  repetition  period  to  calculate 
the  number  of  pulses  received  from  the  target  in  one  scan.  The 
antenna  gain  and  the  number  of  received  pulses  per  scan  are  prime 
parameters  in  the  subsequent  calculations  of  the  signal  power  and 
the  detection  probability. 

The  antenna  equations  as  they  now  exist  can  only  be 
used  for  radars  that  use  the  same  antenna  for  both  transmit  and 
receive.  This,  however,  is  the  most  common  case. 

Signal.  The  signal  power  is  calculated  with  the  standard 
range  equation  using  the  pertinent  parameters  supplied  by  the  input 
cards.  The  signal  power  is  determined  for  the  direct  signal  from 
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the  radar  to  the  target,  the  single  bounce  signal  which  takes  one 
path  to  the  target  via  a  surface  reflection,  and  the  double  bounce 
signal  which  reflects  off  the  surface  to  and  from  the  target.  The 
appropriate  antenna  gain  is  used  in  each  calculation.  That  is, 
in  many  cases  the  angle  taken  by  the  ref lected-signal  path  is 
greater  than  the  half -beamwidth  of  the  antenna  which  puts  it  in 
the  antenna's  sidelobes.  In  this  case  the  signal  power  for  the 
single -bounce  and  double -bounce  signals  are  reduced  from  the 
mainlobe  signal  by  the  given  antenna  sidelobe  ratio  (squared). 

Clutter  and  Receiver-Noise  Power.  The  signal  power 

will  be  divided  by  the  clutter  and  receiver-noise  power  to  form 

the  "signal-to-noise"  ratio  used  in  the  detection-probability 

calculations.  The  receiver  noise  is  calculated  using  the  standard 

equation  involving  Boltzman's  constant,  the  signal  bandwidth,  and 

the  receiver  noise  figure;  rain  backscatter  is  added  to  receiver  noise. 

The  clutter  power,  on  the  other  hand,  involves  numerous 

parameters  such  as  the  antenna  beamwidth  and  the  geometrical 

quantities  shown  in  Figure  3-3.  The  clutter  power  is  calculated 

the  same  way  as  the  target  signal  power,  i.e.,  using  the  radar 

range  equation,  but  using  the  clutter  cross  section  instead  of  the 

target  cross  section.  The  clutter  cross  section  is  equal  to  the 

clutter  area  times  the  normalized  clutter  cross  section  o  ,  which 

o 

was  discussed  in  Chapter  II. 

As  indicated  in  Figure  3-3  the  clutter  area  to  be  used 

in  the  calculation  depends  on  the  range  to  the  target.  The  reason 

'  =* 

for  this  is  that  the  region  on  the  earth's  surface  that  must  be 
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ELEVATION  PLANE 


CLUTTER  INTERCEPTED  BY  PULSE  WIDTH 


fir  TE:  CASE  1  TARGET  COMPETES  WITH  MAINLOBE 
ANDSIDELOBE  CLUTTER 

CASE  2  TARGET  COMPETES  WITH  SIDELOBE 
CLUTTER 


Fig.  3-3  CLUTTER  AREA 
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used  in  calculating  the  clutter  power  is  that  region  which  is 
at  the  same  range  from  the  radar  as  the  target.  When  the  target 
is  low  in  altitude  and  at  a  long  range  then  the  main  beam  intercepts 
both  the  target  and  the  surface  of  the  earth.  In  this  instance, 
called  Case  1  in  Figure  3-3,  the  clutter  area  is  made  up  of  two 
pieces ,  one  in  the  mainlobe  of  the  antenna,  the  other  in  the 
sidelobe .  In  either  case  the  clutter  area  is  bounded  by  the 
pulse  width  in  range. 

When  the  target  is  close  to  the  radar  the  clatter 
region  may  be  at  an  elevation  angle  that  is  out  of  the  main  beam 
of  the  antenna,  such  as  Case  2  in  Figure  3-3.  Even  though  the 
physical  size  of  the  clutter  area  for  this  case  may  be  comparable 
to  that  of  Case  1  the  resultant  clutter  power  received  by  the 
radar  receiver  is  considerably  reduced  due  to  the  low  gain  of  the 
sidelobes . 


Accuracy  Calculations.  The  previous  calculations  can 
be  combined  to  form  the  signal-to-noise  ratio  which  is  important 
in  determining  the  radar's  measurement  accuracy  as  well  as  the 
detection  probability.  Once  the  target  is  detected  the  radar 
must  determine  the  target's  location,  in  range  and  angle,  and  its 
doppler  velocity  or  range  rate.  Since  the  target  echo  exists  in 
the  presence  of  receiver  noise  or  clutter  the  ability  of  the  radar 
to  extract  the  information  is  a  problem  in  the  statistical 
estimation  of  parameters.  The  general  form  of  the  error  for 
measurements  in  range,  range  rate,  or  angle  are  of  the  form3 

aD.  K.  Barton,  Radar  System  Analysis,  (Englewood  Cliffs:  Prentice- 
Hall,  Inc.,  1964),  pp.  38-63. 
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where  E/N  is  equal  to  the  signal-to-noise  energy  ratio  which 
o 

includes  the  observation  time.  The  factor  y  changes  according 
to  whether  range,  doppler  velocity,  or  angle  is  being  measured 
and  according  to  other  factors,  such  as  the  type  of  signal 
modulation.  For  the  range  measurement  the  y  factor  is  basically 
equal  , to  the  bandwidth  of  the  signal;  for  doppler  velocity  it  is 
equal  to  the  length  of  the  signal;  for  angle  measurements  it  is 
equal  to  the  inverse  of  the  beamwidth. 

Thus  for  accurate  measurements  of  the  three  parameters 
it  would  be  desirable  to  have  a  wide  bandwidth  signal,  transmitted 
over  a  long  time  duration,  and  using  an  antenna  with  a  very  narrow 
beamwidth.  In  addition  the  signal-to-noise  ratio  mv3t  be  high. 

The  mathematical  form  of  the  accuracy  equation  is  the 
same  for  both  the  search  accuracies  and  the  tracking  accuracies. 
The  prime  difference  is  that  the  tracking  accuracies  are  usually 
better  due  to  the  longer  time  of  the  measurement.  The  search 
accuracies  are  naturally  lower  then  the  tracking  accuracies  since 
the  time  for  measurement  is  short  due  to  the  finite  time  that  the 
target  remains  in  the  scanning  beam. 

3.4  Detection  Calculations 

The  detection  problem  was  discussed  briefly  in  Chapters 
I  and  II.  It  was  indicated  (See  Figure  2-3)  that  two  detection 
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theories  are  used:  Marcum-Swerling3  and  Shotland4  theory.  The 
Marcum-Swerling  theory  is  a  generalized  theory  for  calculating 
the  probability  of  detecting  a  target,  with  a  fluctuating  cross 
section,  in  the  presence  of  receiver  noise.  It  was  put  into  a 
form,  suitable  for  calculation  by  a  digital  computer,  by  Fehlner6 
and  is  incorporated  in  this  program  as  subroutine  MARCUM.  The 
Shotland  theory  is  an  extension  of  the  Marcum-Swerling  theory  to 
include  the  ffect  of  clutter  as  an  interfering  noise  source  as 
well  as  receiver  noise.  The  Shotland  theory  is  programmed  as 
subroutine  EDDIE. 

Generalized  Detection  Calculations.  The  basic  parameters 
required  to  calculate  the  detection  probability  are  the  false-alarm 
time  from  the  input  cards,  and  the  signal-to-noise  ratio  and  number 
of  received  pulses  per  scan  from  the  radar  calculations.  The 
detection  of  a  target  is  a  decision  process  and  is  based  on 
establishing  a  threshold  level  at  the  output  of  the  receiver.  A 
target  is  assumed  to  be  pre;  ?nt  if  the  signal  out  of  the  receiver 
is  large  enough  to  exceed  th  threshold.  Occasionally  noise  alone 
will  exceed  the  threshold  generating  a  false  alarm.  The  threshold 
level  is  calculated  as  a  function  of  the  noise  power,  the  bandwidth, 

aJ.  I.  Marcum  and  P.  Swerling,  ’’Studies  of  Target  Detection  by 
Pulsed  Radar,"  IRE  Transactions  on  Information  Theory,  vol.  IT-6 
(April,  1960). 

4E.  Shotland,  "False  Alarm  Probabilities  for  Receiver  Noise  and 
Sea  Clutter,"  JHU  Applied  Physics  Laboratory,  BBD-1387 , October  27, 
1964. 

®L,  F.  Fehlner,  "Marcum’s  and  Swerling’s  Data  on  Target  Detection 
by  a  Pulsed  Radar,"  JHU  Applied  Physics  Laboratory,  TG-451, 

July  2,  1962. 
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and  the  number  of  pulses,  to  produce  the  false  alarm  time  specified 
by  the  program  input.  The  probability  of  detecting  a  target  is 
then  the  probability  that  the  target  signal  exceeds  this  pre¬ 
determined  threshold  level  which  varies  up  or  down  according  to 
the  level  of  the  noise  power.  The  procedure  for  calculating  the 

detection  probabilities  involves  evaluating  a  series  of  incomplete 

0 

Gamma  functions,  where  each  function  is  expressed  as  a  series. 

Target  and  Clutter  Statistics.  The  two  detection  theories 
provided  in  this  program  provide  a  fairly  flexible  arrangement  for 
calculating  the  detection  probability  for  targets  and  clutter  with 
various  statistical  characteristics. 

In  the  Marcum-Swerling  theory  incorporated  in  this 
program  five  target  models  may  be  specified.  The  first,  called 
Case  0,  is  a  target  with  a  constant  cross  section  area.  The 
targets  in  the  other  four  cases  have  fluctuating  cross  sections. 

Case  1  and  2  targets  have  a  Rayleigh  amplitude  distribution  and 
apply  to  targets  that  can  be  represented  as  a  number  of  independently 
fluctuating  reflectors  of  about  equal  echoing  area.  Case  1  and  2 
are  often  used  for  aircraft  targets.  Cases  3  and  4  apply  to 
targets  that  can  be  represented  as  one  large  reflector  with  a 
number  of  small  reflectors,  such  as  a  missile.  Cases  1  and  3 
apply  foi  targets  whose  cross  section  fluctuates  slowly,  or  with 
a  period  equal  to  the  scan  time.  Cases  2  and  4  are  for  targets 
whose  cross  section  fluctuates  rapidly,  i.e.  from  pulse-to-pulse . 

3 Ibid. ,  pp.  21-32. 

7Marcum  and  Swerling,  loc .  cit . 
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Since  the  detection  probability  is  higher  for  targets 
which  fall  in  Cases  2  or  4  schemes  are  often  devised  to  insure 
that  the  cross  section  fluctuates  rapidly.  One  way  of  doing 
this  is  to  employ  frequency  jumping  ^ ich  is  accomplished  by 
varying  the  transmitted  frequency  from  pulse-to-pulse  or  within 
a  pulse. 

In  the  Shotland  theory,  as  used  in  this  program, both 
the  target  and  the  clutter  are  of  the  Case  1  variety,  slowly 
fluctuating  with  Rayleigh  amplitude  distributions.  The  clutter 
does  fluctuate  somewhat  faster  than  the  target  depending  on  the 
width  of  the  clutter  spectrum  (calculated  in  the  radar  subroutine) . 
Under  these  conditions,  however,  very  large  signal-to-noise  ratios 
are  required  for  detection.  On  the  other  hand,  if  frequency 
jumping  is  used  with  a  sufficiently  wide  frequency  range  the 
target  changes  to  a.  Case  2  or  4  target  and  in  addition  the  clutter 
spectrum  is  widened  causing  the  clutter  to  fluctuate  randomly 
like  receiver  noise.  Under  these  conditions  the  Mar cum-Swer ling 
theory  can  be  used  which  will  yield  a  much  higher  probability  of 
detection. 

3 . 5  Summary 

The  processing  part  of  the  program  contains  the 
calculations  performed  on  the  input  data  to  ascertain  the 
performance  of  the  proposed  radar  in  the  giyen  situation.  The 
calculations  are  divided  into  three  groups  on  geometry,  radar, 
and  detection. 
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As  indicated  in  the  program  flow  diagram  the  output  for 
a  particular  run  is  stored  in  an  array  and  is  not  printed  until  the 
calculations  for  that  run  are  completed.  This  procedure  has  the 
advantage  of  assuring  that  all  the  results  of  calculations  in 
one  subroutine  are  available  to  other  subroutines.  In  addition, 
it  allows  all  of  the  output  commands  and  formats  to  be  placed  in 
the  main  program. 

The  calculations  and  processing  techniques  discussed  in 
this  Chapter  are  used  to  produce  the  program  output  described  in 
the  next  Chapter . 
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IV.  PROGRAM  OUTPUT 

The  results  of  the  calculations  discussed  in  Chapter  III 
are  provided  as  output  in  either  a  printed  or  digitally  plotted 
format.  The  bulk  of  the  output  is  printed  by  the  main  program, 
as  previously  mentioned,  except  for  the  target  motion  and  cross 
section  decks  which  are  read  and  printed  by  separate  subroutines. 
Two  samples  of  the  program  output  are  included  in  this  Chapter, 
one  on  the  time  simulation  of  the  airborne-radar /missile-target 
example,  proposed  in  Chapter  II,  and  the  other  on  the  range  profile 
for  the  same  radar. 

In  review  of  the  inputs  proposed  in  Chapter  II  the 
airborne-radar /missile-target  example  to  an  airborne  radar  at 
an  altitude  of  30,000  feet  and  at  a  range  of  approximately  43 
nautical  miles  (diagonal  distance  of  30  nautical  miles  on  a 
square)  from  a  ballistic  missile  launched  from  the  surface.  The 
radar  operates  at  a  nominal  frequency  of  5000  megacycles  per 
second  with  a  peak  power  of  500  kilowatts  and  employs  frequency 
jumping  from  pulse-to-pulse .  The  ballistic  missile  flies  in  a 
vertical  plane  for  41  seconds  reaching  an  altitude  of  24,000  feet. 
The  output  of  the  program  in  the  time  simulation  mode  will  cover 
39  seconds  of  the  engagement  and  will  indicate  the  following:  the 
success  or  failure  of  detecting  the  target,  the  length  of  time  from 
missile  launch  required  for  detection,  the  accuracies  of  locating 
the  target,  etc. 
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For  an  example  of  the  range  profile  option  the  target 
is  assumed  to  be  at  a  stationary  altitude  of  30,000  feet  and 
have  a  constant  cross  section  of  .1  square  meters.  The  range 
from  the  target  is  decreased  in  increments  of  2  nautical  miles 
from  a  range  of  78  nautical  miles  with  detection  and  radar 
calculations  being  made  each  increment.  The  single-scan  detection 
probability  and  the  signal-to-noise  ratio  are  plotted  versus 
range  by  the  CALCOMP  Plotter. 

4.1  Time  Simulation  Example 

The  output  for  the  time  simulation  of  the  airborne- 
radar/missile-target  example  is  divided  into  four  parts  for 
discussion  purposes.  The  first  part  is  a  listing  of  the  input 
data.  The  remaining  three  parts  give  the  output  data  arranged 
according  to  computed  constants,  geometry  calculations,  and 
radar  and  detection  calculations . 

4.1.1  Input  Data 

The  data  inputs  are  printed  out  in  three  parts:  target 
motion  deck,  target  cross  section  deck,  and  input  parameter  cards. 
The  target  motion  for  the  ballistic  missile  is  reproduced  in 
Table  4-1.  The  cross  section  of  the  ballistic  missile  is  given 
in  Table  4-2.  The  cross  section  values  are  listed  in  decibels 
above  a  one  square  meter  reference,  for  angles  between  0  and 
36.6  degrees.1  The  values  are  arranged  in  the  format  previously 
described  in  Table  2-2. 

1  Remainder  of  table  to  180°  aspect  angle  not  shown. 
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Table  4-1 

PROGRAM  OUTPUT:  TARGET -MOTION  INPUT 
EXAMPLE  TAKti  dUUST  TRAJECTORY  COORDINATES 


uGWN  RANGE  ALTITUDE 


CARO 

TIME 

POSITION 

SPELU 

ACCELERATION 

PUS  I T I  ON 

SPEED 

ACCELERATION 

1 

0. 

0. 

0. 

0. 

0. 

0. 

25.66 

2 

1  . 

0. 

0. 

0. 

13.0 

25.9 

26.14 

i 

2  • 

0. 

-0. 

0.01 

5/.0 

52.3 

26.59 

4 

3* 

0. 

-0. 

0.02 

11  .0 

79.  1 

27.01 

5 

4. 

0. 

0.1 

0.03 

2  lu.  0 

106.3 

27.40 

6 

3r 

0. 

0.1 

0.03 

330.0 

133.8 

27.74 

7 

6  A 

0. 

0.1 

0.04 

478.0 

161.8 

28.  14 

8 

7  A 

0. 

1.4 

3.49 

654.0 

190.  1 

28.45 

9 

a. 

5.0 

7.3 

8.13 

858.0 

218.6 

28.47 

10 

9  A 

17.0 

16.5 

9.55 

1091.0 

247.2 

28.77 

11 

10. 

38.0 

26.6 

10.73 

1353.0 

276.  1 

29.05 

12 

11. 

70.0 

37.9 

11.84 

1643.0 

305.3 

29.36 

13 

12. 

114.0 

50.2 

12.85 

1963.0 

334.8 

29.73 

14 

13. 

171.0 

63.5 

13.74 

2313.0 

364.7 

30.08 

15 

14. 

242.0 

77.7 

14.53 

2693.0 

395.0 

30.48 

16 

15  . 

327.0 

92.o 

15.23 

3103.0 

425.7 

30.93 

17 

16. 

427.0 

108.1 

15.83 

3545.0 

456.9 

31.39 

18 

17. 

543.0 

124.2 

16.33 

4017.0 

488.5 

31.90 

19 

13. 

675.0 

•  141.1 

17.51 

4522.0 

520.6 

32.27 

20 

19. 

825.0 

159.2 

18.70 

5059.0 

553.0 

32.62 

21 

20. 

994.0 

178.5 

19.89 

5628.0 

585.3 

32.96 

22 

21. 

1183.0 

199.0 

21.11 

6230.0 

618.9 

33.34 

23 

22. 

1392.0 

220.7 

22.31 

6866.0 

652.5 

33.64 

24 

23. 

1624.0 

243.6 

23.52 

7535.0 

686.3 

33.93 

25 

C4. 

1880.0 

267.7 

24.72 

8239.0 

720.3 

34.19 

26 

25. 

2160.0 

293.1 

25.86 

8976.0 

754.  6 

34.27 

27 

26. 

2466.0 

319.5 

26.99 

9748.0 

788.9 

34.33 

28 

27. 

2799.0 

347.0 

28.03 

10554.0 

823.2 

34.18 

29 

28. 

3160.0 

375.6 

29.01 

11394.0 

857.2 

33.92 

30 

29. 

3550.0 

405.0 

29.84 

12268.0 

890.9 

33.33 

31 

30. 

3970.0 

435.1 

30.26 

13176.0 

923.6 

31.94 

32 

31. 

4420.0 

465.4 

30.12 

14115.0 

954.4 

29.49 

33 

32. 

4900.0 

495.8 

30.89 

15083.0 

983.5 

28.99 

34 

33. 

5411.0 

527.2 

31.85 

16081.0 

1012.4 

28.69 

35 

34. 

5954. C 

559.7 

33.13 

17108.0 

1041.4 

29.37 

36 

33. 

6530.0 

593.6 

34.54 

18165.0 

1071.1 

30.07 

37 

36. 

7141.0 

628.9 

36.03 

19251.0 

1101.6 

30.87 

38 

37. 

7787.0 

665.7 

37.61 

20368.0 

1132.9 

3 1 »  75 

39 

38. 

8471.0 

704.2 

39.28 

21517.0 

1165.1 

32.67 

40 

39. 

9194.0 

744.3 

40.87 

22698.0 

1198.3 

33.57 

41 

40. 

9959.0 

786.0 

42.60 

23914.0 

1232.3 

34.  58 
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The  input  parameters  that  are  used  to  describe  the 
initial  radar-target  geometry,  the  clutter  and  noise  conditions, 
the  radar  characteristics,  and  the  output  options,  are  reproduced 
in  Table  4-3.  These  parameters,  which  were  described  in  Chapter 
II,  are  given  the  values  listed  in  the  nominal-value  column.  The 
nominal  vlluep  are  multiplied  by  the  number  given  in  the  conversion 
factor  column  to  agree  with  the  units  in  the  definition  column. 

4,1,2  Computed  Constants 

The  first  set  of  output  parameters  are  constant  with 
respect  to  time  and  are  therefore  calculated  only  once  during 
each  run.  The  results  of  these  calculations  are  listed  with 

tions  in  Table (4-4.  Note  that  the  horizontal  scan  rate  of 
the  antenna  was  computed  as  800  degrees  per  second.  This  implies 
that  the  airborne  radar  must  employ  an  electronically  scanned 

antenna  since  mechanical  scanning  cannot  be  performed  at  an  800 
degree/second  rate. 


4,1,3  Geometry  Output 

The  results  of  the  geometry  calculations  occur  on 
three  pages  of  the  printed  output.  The  first  page  is  for  the 
ranges  and  range  rates  between  the  radar,  target,  reflection  point 
and  clutter  spot,-  the  second  page  is  for  the  angle  calculations; 
and  the  third  page  is  for  the  cross-section  calculations. 

The  definitions  of  all  the  geometry  output  parameters 


are  listed  in  Table  f-5  in  two  parts: 
angle  calculations  and  part  II  for  the 


part  I  for  the  range  and 
cross-section  calculations. 
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Name 

LAMBDA 

AR 

GAIN 

BETAH 

BETAV 

FSCH 

TAU 

N 

PO 

PAVG 

PNOISE 

PFA 


Table  4-4 

PROGRAM  OUTPUT:  COMPUTED  CONSTANTS  FOR 
AIRBORNE -RADAR/M I S SI LE -TARGET  EXAMPLE 


Definition 
Wavelength  in  feet 

Effective  area  of  antenna  in  square  feet 
Main lobe  gain  in  dB 
Horizontal  beamdwidth  in  degrees 
Vertica*  Deamwidth  in  degrees 
Horizontal  scan  rate  in  degrees/sec 
Compressed  pulse  length  in  seconds 

Number  of  hits/scan 

\ 

Power  received  at  initial  range  in 
watts/square  feet 

Average  transmit  power  in  watts 

Noise  power  in  dB  above  one  milliwatt 
plus  rain  backscatter  power  when  RAIN  = 

False  alarm  probability 


Computed 
1  Value 

.  1967 

20.66 

38 

2.0 

2.0 

800 

1.0  x  10~7 

5 

2.0  x  10"12 
5  x  103 

-99 

1.4  x  10'8 


1. 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 
Silver  Sprino  Maryland 


RDOT 

R1D0T 

R2D0T 

TDIR 

DELT 


Radar- to-target  range  in  feet 
Radar-to-ref lection-point  range  in  feet 
j Target-to-ref lection-point  range  in  feet 


Respective  range  rates  in  feet/sec. 


Propagation  time  of  direct  signal  in  secs. 

Time  difference  between  one-way  direct 
signal  and  one-way  reflection  signal  in  secs 


Page  2:  Geometry  Ancle  Calculations 


ELAAT  Elevation  angle  of  radar/target  line  in 

degrees 

ELAAB  Elevation  of  radar/reflection-point  line 

in  degrees 

ELAAC  Elevation  of  radar /clutter-point  line  in 

degrees 

AZAAT  Azimuth  angle  of  radar/target  line  in 

degrees 

ALPHAB  Angle  of  incidence  of  radar/reflection-point 

line  with  respect  to  tangent  at  reflection 
point,  in  degrees 

ALPHAC  Angle  of  incidence  with  respect  to  tangent 
at  clutter  point,  in  degrees 

THETA (1)  Angle  .between  target  plane  of  rotation  and 

radar  polarization  plane  in  degrees 
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Table  4-5 

DEFINITIONS  OF  GEOMETRY  OUTPUT  PARAMETERS,  PART  II 


Page  3:  Geometry  Cross-Sections  Calculations 

ALP(l)  Direct-path  aspect  angle  in  degrees 

SIG(1,1)  Target  cross  section  for  direct  path  and 

horizontal  polarization,  in  square  feet  (POL 

SIG(1,2)  Target  cross  section  for  direct  signal 

and  vertical  polarization,  in  square  feet 

ALP(2)  Single-bounce-path  aspect  angle  in  degrees 

SIG(2 , 1)  Target  cross  section  for  single-bounce  path 

and  horizontal  polarization,  in  square  feet 

SIG(2,2)  Target  cross  section  for  single-bounce  path 
and  vertical  polarization,  in  square  feet 

ALP (3)  Double-bounce  path  aspect  angle  in  degrees 

SIG(3,1)  Target  cross  section  for  double-bounce  path 
and  horizontal  polarization,  in  square  feet 

SIG(3 , 2)  Target  cross  section  for  double-bounce  path 
and  vertical  polarization,  in  square  feet 
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The  geometrical  qua<ntitites  that  these  definitions  describe  were 
previously  illustrated  in  Figure  3-2. 

The  actual  output  for  the  geometry  range  calculations 
for  the  airborne-radar/missile-target  example  is  reproduced  in 
Table  4-6.  This  Table  indicates  that  the  range  R  changes  only 
16,000  feet  (2.7  nautical  miles)  from  the  initial  range,  increasing 
from  zero  to  36  seconds  and  decreasing  thereafter. 

The  geometry  angle  calculations  are  shown  in  Table  4-7. 
Note  that  all  of  the  elevation  angles  are  negative  indicating 
that  the  antenna  beam  is  pointing  down  throughout  the  engagement. 

The  last  page  of  the  geometry  output  showing  the  target 
cross-section  values  for  the  different  signal  paths  is  reproduced 
in  Table  4-8 . 

4.1.4  Radar  and  Detection  Output 

The  results  of  the  radar  and  detection  calculations  are 
presented  on  two  pages.  The  first  page  includes  the  output 
parameters  associated  with  the  detection  probabilities  and  the 
accuracy  calculations  which  are  defined  in  Table  4-9.  The 
second  page  presents  the  output  parameters  associated  with  the 
signal  and  clutter  power  level  calculations  which  are  defined  in 
Table  4-10.  The  corresponding  two  pages  of  output  for  the  airborne- 
radar/missile-target  example  are  shown  in  Tables  4-11  and  4-12. 

Some  of  the  most  important  conclusions  revealed  by  these  two 
tables  are  indicated  below. 

The  first  measurable  single  scan  de\  tion  probability 
PD1,  equal  to  .403,  occurs  11  seconds  after  the  missile  launch. 
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Table  4-8 

PROGRAM  OUTPUT:  GEOMETRY  CROSS-SECTION  CALCULATIONS 

S I  G(  1*2)  ALP  ( 2 1  SIG(  2»  1  )  SIG(2,2)  ALP<3i  SIG(3,1)  SIG 
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PD1 

PDCUMl 

PDT 


PCUMT 


CASE 

VB 

S/N 

C/N 

MC 

FD1 

DOPDIF 

TDIF 

DELR 

DELD 

DELA 

SIGNV 

SIGNH 

SIGT 

ZT 


Table  4-9 


DEFINITIONS  OF  DETECTION  PROBABILITIES  AND 
RADAR  ACCURACY  CALCULATIONS 


Single-scan  detection  probability  for  the  direct 
path  signal 

Cumulative  detection  probability 

Detection  probability  for  direct  and  reflected 
paths  (when  coincide) 

Cumulative  detection  probability  of  combined 
signal 

Marcum-Swerling  f luctuating-target  model  (0-4) 

Threshold  level  above  receiver  noise 

Signal-to- (receiver)  noise  ratio  in  dB 

Clutter-to- (receiver)  noise  ratio  in  dB 

Number  of  correlated  clutter  pulses 

Absolute  target  doppler  frequency  in  cps 

Difference  in  doppler  frequencies  between  direct 
and  reflected  paths,  in  cps 

Time  difference  between  direct  and  reflected 
signals  in  microseconds 

Search  range  accuracy  in  feet 

Search  doppler  accuracy  in  cps 

Search  angular  accuracy  in  milliradians 

Angular  tracking  accuracy  in  vertical  plane 

Angular  tracking  accuracy  in  horizontal  plane 

Total  angular  tracking  accuracy  in  milliradians 

Altitude  of  target  in  thousands  of  feet 


-  61 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 
Silver  Spring.  Maryland 


Table  4-10 

DEFINITIONS  OF  SIGNAL  AND  CLUTTER 
POWER  LEVEL  CALCULATIONS 


SDD 

Signal  power  for  direct  path  in  dBm 
(decibels  above  one  milliwatt) 

SDB 

Signal  power  for  single-bounce  path  in 

dBm 

SBB 

Signal  power  for  double-bounce  path  in 

dBm 

CLOM 

Mainlobe  clutter  power  in  dBm  (surface 

clutter) 

CLOS 

Sidelobe  clutter  power  in  dBm 

CLO 

Total  resultant  clutter  power  in  dBm 

SIGO 

Normalized  clutter  cross  section  in  dB 

PHA 

Phase  difference  between  direct  and  bounce 
path  in  cycles 

PHADEL 

Phase  change  during  one  scan  in  cycles 

REFL 

Reflection  coefficient 

S/(C+N) 

Sighal-to-clutter  plus  noise  ratio  in  dB 
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Reliable  detection  of  the  missile,  however,  does  not  occur  until 
19  seconds  when  the  detection  probability  equals  unity.  From  19 
seconds  until  the  end  of  the  engagement  the  detection  probability 
fluctuates  between  unity  and  .03.  The  latter  probability  value 
occurs  at  36  seconds  where  the  missile  cross  section  drops  to  its 

smallest  value  of  .42  square  feet. 

The  missile  is  not  detected  reliably  until  19  seconds, 
because  of  the  high  clutter  conditions .  As  shown  in  Table  4—12 
the  clutter  power  level  CLO  is  composed  of  clutter  echoes  entering 
the  mainlobe  of  the  antenna  until  18  seconds.  Aft*  ’  18  seconds 
the  clutter  level  decreases  abruptly  because  at  this  point  the 
clutter  echoes  enter  via  the  antenna  sidelobes  only .  Thus  the 
airborne  radar  in  this  particular  engagement  cannot  detect  the 
missile  in  a  background  of  mainlobe  clutter,  but  must  wait  until 
the  missile  is  high  enough  for  reliable  detection  in  sidelobe 
clutter . 

The  effect  of  the  high  mainlobe  clutter  return  is  also 
reflected  in  the  signal-to-noise-plus-clutter  ratio,  S/(C+N)  in 
Table  4-12.  As  shown  the  ratio  is  less  than  unity  below  19 
seconds  and  greater  than  unity  with  fairly  high  values  above  19 
seconds . 

After  the  radar  detects  the  missile  at  19  seconds  its 
location  is  determined  in  the  search  mode  to  an  accuracy  of  two 
feet  in  range,  358  KHz  in  doppler,  and  four  milliradians  in  angle. 
The  poor  doppler  accuracy  would  force  this  particular  radar  to 
determine  velocity  by  measuring  the  range  rate.  The  angle  accuracy 
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of  the  search  radar  in  this  case  is  high  enough  for  the  radar  to 
begin  tracking  with  an  angular  accuracy  equal  to  one  milliradian 

4.2  Range  Profile  Example 

The  time  simulation  mode  analyzed  the  performance  of 
the  airborne  radar  in  detail  in  detecting  and  locating  the  missile, 
launched  approximately  43  nautical  miles  from  the  radar.  The 
range  profile  output  discussed  here  will  indicate  the  detection 
performance  of  the  same  airborne  radar  but  this  time  against  a 
simplified,  generalized  target. 

The  target  for  the  range  profile  run  is  described  by 
the  input  parameters  listed  in  Table  4-13.  As  shown  the  target 
analyzed  will  remain  at  an  altitude  of  30,000  feet  and  will  have 
a  constant  cross  section  of  .1  square  meters  (-10  dB/square  meters), 
The  three  input  parameters  shown  are  the  only  ones  that  are  changed 
from  the  previous  time-simulation  run.  They  specify  calculations 
to  be  performed  every  two  nautical  miles  from  zero  to  a  maximum  of 
78  nautical  miles. 

The  results  are  shown  in  Tables  4-14  and  4-15  in  the 
same  format  used  in  the  time  simulation  output  except  that  the 
target  range  is  incremented  instead  of  the  time  variable.  The 
radar  itself  is  operating  in  the  same  mode  of  operation  as  used 
in  time  simulation  run,  that  is,  searching  for  the  target  over 
a  26  by  30  degree  sector  in  one  second.  This  produces  five  return 
pulses  which  are  integrated  by  the  radar  to  produce  the  detection 
probabilities  listed  in  Table  4-14.  The  corresponding  signal-to- 
noise-plus-clutter  ratio  is  listed  versus  range  in  Table  4-15. 
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The  detection  probability  PD1  and  signal-to-noise-plus- 
clutter  ratio  are  generally  the  most  important  quantities  in  the 
range  profile  ootion  and  are  therefore  plotted  with  the  CALCOMP 
plotter  as  shown  in  Figure  4-1  End  4-2.  These  figures  indicate 
that  the  nJrborne  radar  would  reliably  detect  a  .1  square  meter 
target  at  an  altitude  of  30,000  feet  out  to  a  range  of  about  40 
nautical  miles.  The  target  altitude  of  30,000  feet  is  significant 
in  that  the  clutter  returns  enter  the  receiver  via  the  antenna 
sidelobes.  If  the  target  was  much  lower  in  altitude  the  clutter 
would  be  in  the  mainlobe  of  the  antenna  which  would  preclude 
target  detection. 
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SINGLE-SCAN  DETECTION  PROBABILITY 


SIGN AL-TO-Cl UTTER. PLUS-NOISE  RATIO  (decibel*) 
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V.  SUMMARY  AND  CONCLUSIONS 

The  digital  computer  program  for  radar  analysis  and 
simulation  has  been  described.  In  this  chapter  some  of  the  high¬ 
lights  will  be  summarized  and  some  areas  for  future  development 
will  be  discussed. 

5 . 1  Summary 

The  report  started  with  a  description  of  the  basic 
radar  problem  to  be  analyzed  in  this  study:  the  detection, 
acquisition  and  tracking  performance  of  a  generalized  radar 
against  a  single  moving  target  in  a  clutter  environment. 

Chapters  II  through  IV  described  the  program  inputs,  the  process¬ 
ing,  and  the  program  outputs,  respectively.  Chapter  II  on  inputs 
described  the  manner  in  which  the  radar  problem  is  defined  by  the 
user.  A  flexible  but  detailed  arrangement  was  indicated  for 
describing  target  motion,  target  cross  section,  initial  radar- 
to-target  geometry,  radar  characteristics,  clutter  and  noise 
environment,  and  program  options  Any  of  the  input  parameters 
may  be  changed  for  performing  multiple  runs  on  any  one  job. 

The  processing  or  calculations  performed  by  the  program, 
as  described  in  Chapter  III,  were  shown  to  be  arranged  as  a  time 
simulation  of  the  radar-target  engagement.  The  geometry  calcu¬ 
lations  are  performed  with  vector  algebra  permitting  an  accurate 
simulation  of  the  radar  and  target  motion,  and  calculation  of  the 
pertinent  geometrical  quantities  between  the  radar,  clutter,  and 
reflection  points.  The  radar  calculations  are  performed  each 
time  step  yielding  signal-to-noise  ratios,  clutter  characteristics 
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and  search  and  tracking  accuracies.  The  detection  probabilities 
are  calculated  for  many  target  and  noise  models,  and  include  the 
cumulative  detection  probabilities.  A  program  option  permits  the 
target  model  to  be  greatly  simplified  for  performing  radar 
comparisons  with  the  target  motion  and  cross  section  held 
constant. 

A  few  samples  of  the  program's  output  were  presented 
in  Chapter  IV.  In  the  time  simulation  mode  the  output  is  arranged 
as  a  function  of  time  with  calculations  being  made  every  time  the 
radar  antenna  scans  past  the  target.  The  output  parameters  are 
quite  extensive  and  include:  single-scan  and  cumulative  detection 
probabilities  for  the  direct  and  surface-reflected  signals;  search 
accuracies  in  range,  doppler,  and  angle;  angular  tracking  accuracies 
clutter  statistics  and  power  levels;  cross  section  variations; 
s ignal-to-noise  ratios;  and  many  geometrical  parameters. 

In  the  range  profile  mode  the  target's  motion  and 
cross  section  'are  held  invariant  to  permit  comparisons  of  radar 
performance o  The  target  model  is  grossly  oversimplified  in  this 
mode  which  is  necessary  to  analyze  the  performance  of  a  radar, 
or  radars,  for  a  broad  range  of  target  cross  sections,  target 
statistics,  clutter  conditions,  and  distances  from  the  radar. 

That  is,  the  boundaries  of  the  radar's  performance  can  be  found 
by  operating  in  the  range  profile  mode  with  calculations  being 
made  for  all  probable  cross  section  values,  clutter  conditions, 
etc.,  without  regard  to  one  particular  radar-target  engagement. 
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Thus,  the  range-profile  mode  is  used  for  calculating 
general  radar  performance,  while  the  time-simulation  mode  is  for 
detailed  radar  performance. 

The  actual  program  was  written  in  FORTRAN  II  in  a  manner 
to  achieve  considerable  flexibility  and  growth  capability.  This 
was  accomplished  by  dividing  the  program  into  22  subroutines 
controlled  by  one  main  program.  Thus,  if  any  changes  are  required, 
or  desired,  they  may  be  inserted  into  the  subroutine  in  which  they 
aPPly  without  affecting  the  status  of  the  other  subroutines. 

5.2  Areas  for  Future  Development 

The  program  as  it  now  exists  is  a  fairly  complete 
representation  of  the  radar-target  engagement.  There  are  some 
areas,  however,  which  deserve  future  attention;  some  of  these 
were  indicated  in  previous  chapters. 

Antenna.  There  are  two  antenna  characteristics  that 
are  simplified  in  the  present  program  that  could  be  changed  to 
achieve  a  more  accurate  representation  of  an  actual  antenna.  They 
are  the  antenna  beam  pattern  and  the  antenna  beam  motion. 

Presently,  the  antenna  beam  pattern  has  a  constant 
mainlobe  gain  ever  the  solid  angle  defined  by  the  half  power 
points,  and  a  constant  sidelobe  gain  for  all  other  angles.  In 
the  future,  it  may  be  desirable  to  modify  the  radar-calculation 
subroutine,  or  add  an  additional  subroutine,  to  incorporate 
antenna  gain  variations  with  angle. 

The  motion  of  the  antenna  beam  is  presently  arranged 
to  scan  the  given  angular  sector  in  the  allotted  scan  time  with 
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the  target  assumed  to  be  in  the  center  of  the  sector.  A  more 
accurate  representation  of  the  antenna  beam  motion  can  be  incor¬ 
porated  by  having  the  antenna  beam  follow  a  predetermined  scan 
pattern  starting  from  an  initial  position  specified  on  the  input 
or  chosen  randomly  by  Monte  Carlo  techniques.  This  method  would 
change  the  time  increment  of  the  output.  That  is,  the  calculations 
of  the  radar  detection  probabilities  would  be  performed  only  at 
those  times  in  which  the  target  is  in  the  antenna  beam.  This  may, 
or  may  not,  occur  once  each  sector-scan  time  as  in  the  current 
program. 


Acquisition  Mode.  Presently  the  transition  from  the 
search  mode  to  the  tracking  mode  is  assumed  to  occur  instantaneously. 
This  somewhat  unrealistic  assumption  can  be  corrected  by  calculating 
the  time  required  to  point  the  antenna  at  the  target.  The  ca.-i- 
lation  can  be  programmed  using  techniques  by  Barton1  and  quantities 
available  in  the  present  program;  namely,  the  time  required  to 
scan  the  search  sector,  the  ratio  between  the  size  of  the  search 
sector  and  the  tracking  beam,  and  the  signal-to-noise  ratio. 


Jamming.  The  present  program  does  not  consider  jamming. 
However,  radar  performance  in  the  presence  of  certain  types  of 
Jamming,  such  as  wide-band  noise  jammers,  could  be  added  in  a 
straightforward  manner  by  modifying  the  signal-to-clutter-plus- 
noise  ratio  to  include  the  noise  power  received  from  a  jammer. 

XD.  K.  Barton,  Radar  System  Analysis.  Prentice-Hall ,  Inc., 
Englewood  Cliffs,  New  Jersey,  1964,  Section  14.2, 
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Other  Areas.  There  are  evidently  other  areas  of  the 
program  that  could  be  improved  or  other  simulations  that  could 
be  added,  to  represent  more  accurately  the  real  situation.  Some 
of  these  areas  are  an  improved  simulation  of  the  tracking  mode, 
clutter  model  for  land  (sea  clutter  model  and  rain  clutter  already 
included) ,  multiple  target  detection,  sequential  detection,  and 
target  motion  not  confined  to  a  vertical  plane. 
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APPENDIX  A 

LITERATURE  REVIEW 

The  literature  was  searched  for  digital  computer 
programs  and  radar  simulation  techniques  applicable  to  this 
Study.  The  pertinent  literature  can  be  conveniently  grouped 
into  the  following  categories:  war  game  simulations,  computer 
programs  for  radar  analysis,  and  detection-probability  algorithms. 

War  Game  Simulations.  Numerous  programs  exist  which 
simulate  the  radar  only  to  the  extent  of  determining  the  outcome 
of  a  war  game .  Examples  of  war  game  simulations  involving  radars 
are  programs  written  by  Andrus  1  ,  and  Kennard8 .  The  program  by 
Andrus  is  designed  to  simulate  the  interactions  between  surface- 

i 

to-air  missile  systems  and  aircraft.  The  different  events  of 
the  game  are  assigned  probability  values  and  during  the  operation 
of  the  program  the  occurrence  of  an  event  is  determined  by  comparing 
the  given  probability  with  an  internally  generated  random  number. 

The  gross  parameters  of  the  system  are  fed  into  the  program  along 
with  the  number  and  positions  of  the  aircraft.  The  program 
determines  missile  firing  doctrines,  optimum  number  and  placement 
of  the  radars,  and  target  kill  probabilities. 

XA.  F.  Andrus,  ”A  Computer  Simulation  for  the  Evaluation  of 
Surface-to-Air  Missile  Systems  in  a  Clear  Environment,” 

Naval  Postgraduate  School,  TR/RP  No.  67,  June  1966. 

8P.  H.  Kennard,  "Effectiveness  Studies  of  Manned  Airborne  Weapon 
Systems  by  Digital  Simulation”.  Abstract:  Instruments  and 
Control  System,  vol.  32  (July,  1961),  pp.  2171-1272. 
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These  types  of  war  game  simulations  are  considerably 
different  from  the  simulation  in  the  Radar  Analysis  Program.  The 
only  radar  parameters  involved  in  the  Andrus3  program  are  the 
search  radar  maximum  range,  the  probability  of  detecting  a  target, 
and  the  time  required  to  acquire  a  target  with  the  tracking  radar. 
These  are  performance  parameters  of  a  radar  and  in  the  Radar 
Analysis  Program  they  are  calculated,  as  opposed  to  being  provided 
as  inputs  in  war  game  simulations. 

Computer  Programs  for  Radar  Analysis.  Three  programs 
were  found.  They  will  be  discussed  in  order  of  their  increasing 
similarity  to  the  Radar  Analysis  Program. 

The  first  program  by  White  and  James4  is  an  example  of 
how  a  digital  computer  program  can  be  used  to  simulate  or  model 
a  radio  communication  system,  which  is  in  man y  ways  similar  to  a. 
radar  system.  Their  specific  problem  was  to  model  the  transmitter, 
receiver ,  and  atmospheric  phenomena  in  order  to  study  interference 
problems.  This  involved  deriving  mathematical  equations  to  describe 
the  operations  of  each  significant  parameter  of  the  transmitter, 
receiver,  and  propagation  phenomena.  Some  of  the  techniques  which 
are  used  for  modeling  the  antenna  gain  pattern  and  the  propagation 
phenomena  are  the  same  as  those  used ' for  the  radar  case  and  can 
be  incorporated  into  the  Radar  Analysis  Program  if  desired. 

aAndrus ,  ojd.  cit .  ,  p.  4. 

D.  R.  J.  White  and  W.  G.  James,  "Digital  Computer  Simulation 
for  Prediction  and  Analysis  of  Electromagnetic  Interference," 

IRE  Transactions  on  Communications  Systems,  vol.  CS-9  no  2 
(June,  1961),  pp.  148-159. 
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In  the  second  example  of  a  program  for  radar  analysis, 

Wan  briefly  describes  a  completely  analytical  approach  to  a 
radar  program  that  he  says  will  be  developed.6  Wan's  technique 
is  to  isolate  parts  of  the  problem,  such  as,  the  transmitter, 
receiver,  antenna,  noise  power,  weather,  and  target  and  to  write 
a  mathematical  model  of  each  in  signal-space  notation.  The 
transmitter,  for  example,  can  be  modeled  strictly  in  the  signal 
or  frequency  domain  by  representing  it  by  the  Fourier  transform 
of  the  transmitted  time  waveform.  The  model  for  the  target,  on 
the  other  hand,  is  described  in  both  signal  and  space  notation; 
the  cross  section  as  a  function  of  the  transmitted  frequency 
describes  the  target  in  the  signal  (or  frequency)  domain,  while 
the  range  of  the  target  from  the  radar  describes  the  target  in 
the  space  domain.  These  mathematical  models,  or  transfer  functions, 
of  the  different  parts  of  the  problem  are  combined  to  simulate  a 
particular  radar-target  engagement.  For  example,  if  the  radar  is 
to  operate  in  rain  the  transfer  function  which  synthesizes  the 
rain  would  be  included. 

The  last  program  discussed  in  this  section  has  the  most 
similarity  to  the  Radar  Analysis  Program.  This  program,  by  Boothe, 
is  used  to  determine  the  performance  of  an  acquisition  radar  through 
application  of  radar  detection  probability  theory.8  Both  programs 


L.A.  Wan,  "Tactical  System  Radar  Signal-Space  Model,"  Ninth 
Conference  on  Military  Electronics.  September  22-24,  1965, 
sponsored  by  IEEE,  pp.  13-17. 


R.  R.  Boothe,  "A  Digital  Computer  Program  for  Determining  the 
Performance  of  an  Acquisition  Radar  through  Application  of  Radar 
Detection  Probability  Theory,"  Advanced  Systems  Laboratory,  U.  S. 
Army  Missile  Command,  No.  RjpjpTR-64-2 ,  December  31,  1964. 
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are  arranged  as  time  simulations  of  a  single  radar,  single 
target  engagement.  In  addition  both  programs  simulate  the  target 
motion  and  cross  section  by  means  of  external  input  for  maximum 
flexibility  and  use  the  standard  radar  equation  for  calculating 
signal-to-noise  ratio.  Some  of  the  notable  differences  between 
the  Boothe  program  and  the  Radar  Analysis  Program  are  in  the 
calculation  of  the  noise  power  from  all  sources  and  in  the 
calculation  of  the  detection  probability.  The  Boothe  program 
ignores  clutter  but  does  have  an  option  for  including  jammer 
noise.  The  Boothe  program  incorporates  a  detection  theory 
developed  by  Marcum7  in  1947  which  is  a  method  for  calculating 
the  detection  probability  for  a  constant  cross  section  target 
in  a  background  of  receiver  noise.  In  1960  Marcum  in  conjunction 
with  Swerling  extended  the  theory  to  include  fluctuating  targets.8 
This  theory  is  used  in  the  Radar  Analysis  Program  in  a  computer 
program  developed  by  Fehlner.® 

Target  Detection  Programs.  Two  other  sources  were  found 
which  contain  digital  computer  programs  for  calculating  detection 
probability.  Both  of  these  programs  use  Marcum  and  Swerling 
detection  theory. 


7J.  I.  Marcum,  "A  Statistical  Theory  of  Target  Detection  by  Pulsed 
Radar,"  The  Ra.nd  Corporation,  RM-754,  December  1,  1947. 

8J.  I.  Marcum  and  P.  Swerling,  "Studies  of  Target  Detection  by 
Pulsed  Radar,"  Institute  of  Radio  Engineers  Transactions  on 
Information  Theory.  Vol.  IT-6.  (April  1960). 

SL.  F.  Fehlner,  "Marcum's  and  Swerling' s  Data  on  Target  Detection 
by  Pulsed  Radar,"  The  Johns  Hopkins  University  Applied  Physics 
Laboratory,  TG-451,  July  2,  1962. 
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The  first  program,  by  Kirkwood,  calculates  the  cumulative 
detection  probability  for  targets  approaching  a  scanning  radar.10 
In  Kirkwood’s  program  only  two  out  of  the  four  Marcum-Swerling 
target  models  are  used,  and  the  single  scan  detection  probabilities 
are  not  available.  The  program  is  available  in  the  reference  and 
is  written  in  FORTRAN. 

The  other  computer  program  for  calculating  detection 
probability,  developed  by  Nolen,  et  al.11  is  very  similar  to  the 
work  performed  by  Fehlner.  Both  of  these  programs  were  produced 
for  the  prime  purpose  of  generating  a  set  of  general  charts  to 
be  used  by  radar  analysts  for  calculating  the  detection  probability 
manually.  The  programs  were  run  for  a  large  number  of  cases  to 
cover  most  conditions.  The  charts  produced  by  the  Nolen  program, 
however,  are  only  available  for  20  integrated  pulses  while  Fehlner's 
charts  are  available  for  up  to  3000  integrated  pulses.  Th.»  program 
developed  by  Fehlner  was  arranged  in  a  subroutine  and  is  the  main 
detection  probability  algorithm  in  the  Radar  Analysis  Program. 


l0P.  K.  Kirkwood,  ’’Radar  Cumulative  Detection  Probabilities  for 
Radial  and  Nonradian  Target  Approaches,"  The  Rand  Corporation, 
RM  4653-PR,  September  1965. 


11 J.  C.  Nolen,  et  al.  /'Statistics  of  Radar  Detection,"  Th^  Bendix 
Corporation,  Baltimore,  Maryland,  February  1966. 
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Table  B-l 

SUBROUTINE  NAMES  IN  PROGRAM 
LISTING 


Section  No. 

1.0 

1.1 

1.2 

2.0 

2.1 

2.2 

2.3 

2.4 

2.5 

2.6 

2.7 

2.8 

2.9 


DAVE 

3.0 

EDDIE 

4.0 

PFA(f) 

4.1 

PFAC(f) 

4.2 

MARCUM 

4.3 

DGAM(f) 

4.3.1 

DEVAL(f) 

4.3.2 

GAM(f) 

4.3.3 

EVAL(f) 

4.3.4 

SUMLOG(f) 

4.3.5 

MAINP 

TARGIN 

CROSIN 

GEOM 

TARGET 

AIRCFT 

UNIT 

CROSS 

TRIAD 

MULT 

DOT 

ANGLE 

RATSCT 


1.0 

CRAU Ak 


I  KAuAk  analysis  program 

OIM.LNbluN  0(0) 

OIMlNSIuN  GPK(iU)  -  ' 

OlMtNSiUN  XX  A  a  (  A  ) ,  YXAX(2),  XYAX(2),  YY^T2i 

01  MLkS  ION  XPH{  <£i>  »  YPH  ( 2)  ,  XPVI2),  YP  VJ/tn 

01  ME Nil UN  EPH(  322,0) ,CkUSEC(  18  01  .  2)yGI  (  30 , 32 )  *H(  80 , 4u 

CgM.MiLN  1  I  ML,  KSlcf-  j  KHT  ,  LPH,  CRCSEC  » H  »CbM 

WO  IK  =  6u 

Uu  b  1=1,100 

CoM( 1 )  =0. 


Jf 

fh X(2),  YY^fti 
2),  YPV^fT 


KtNC  =  1 

LnCPLT  =  (  +6HtN0PL  I  I 
KHt  0  =  0 


CGM( 100) 


CALL  1  ttKbl  N(  TMAX) 

CALL  CKUdlN 

0  k  t  AD  I  NP  U  J  JAPt  b*i2»NDUMA*  NDUMb  ,KC,SYM,A,B,(L)(  1)  ,  1  =  1,6) 

12  1-ukMAI  (  1  i  ,  1  2 ,  1 3  ,  At> ,  2 1 1  b  •  d  ,  o A6  ) 
lF(NDJMA)  iul » 102, 101 

10 1  Call  Takgin(Tmax) 

10 2  lF(NDUMD)  lu3»  lu4 »  ic3 

103  CALL  CRUSIN 
Gu  TC  9 

luA  IF(NUUMa)  9,lob»S 

105  1 F ( KC)  ib,lb,loo 

106  1 F (  kHE  0 )  13 , lu7 , 13 

107  rtKlTt  DO  TP  UT  TAPE  6,10 

10  FORMAT  tlHl  bnX  cLHALblS  INPUT  PARAMETERS  //18H  II  12  13  SYMBOL 

1  7X  13HNUMINAL  VALUE  3X  1 7HCONV6RS 1  ON  FACTOR  8X 

2  luHOEF INI T ION  ///) 

KHfcL)  =  1 

12  wklTt  UUlPUT  TAPE  o, 14,N0UMA,N0UMb,KC,SYM,A,B, (D( I ) , 1=1,6) 

14  FORMAT  <2X,Ii,I3,I4,2X,A6,2E2C.8,2X,6Ao) 

CuM(KC)  =  A*B 

GO  TO  9 

15  CONTINUE 

IF tSYM-ENOPLT)  17,16,17 

16  KENC  =  2 

17  KHfcU  =  0 
KHI  =  0 

CuM ( 31 )  =  CUM(2)*CUSP<CUM<3))*C0SF(C0M<4)) 

CUM ( 32 )  =  CUM{ 2 ) *COSF ( CGM{ 3 ) ) *  SI NF ( C0M(4 ) } 

COM ( 33 i  =  CUM( 2) *S INF (COM I  3 ) ) 

COM  (  34  )  =  SI  INF  (COM  I  lu)  ) 

CUM  I  3b )  =  COSF ( COM ( 10) ) 

I F (COM ( 6 1 ) )  28,26,18 

18  TIME  =  TMAX 
SNAML  =  ( +5HRANGE ) 

TNAMt  =  (+4HJN  Ml) 

NHT  =  <CUM(62)-C0M(5))/C0M(61)  +  1.01 
NkEM  =  NHT 

IS  IF (NREM-NOIM )  21,21,22 

21  NSTEP  =  NREM 
NkEM  =  0 

GU  TO  23 

22  NSTEP  =  NDIM 
NkEM  =  NREM-NOIM 

23  Ou  20  K=l, NSTEP 
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CUHtil  =  CUK(ol)*FLUATF(KHl)  +  COM%  0 ) 

KbIEP  =  K 
KHT  =  KhT  +  1 
CALL  GEuM 
CALL  GAVt 
CALL  EbolE 

20  H( K , 1 )  *  COi'l  (o3)  /oo80. 

NL  *  4 

GU  TL  iUi 

2b  CGM63)  =  CUMi 5 ) 

SNAML  =  (+4HTIME) 

TNAMt  =  (♦3HSEC) 

HTSS  =  CUM29) 

NHT  =  MlNIFi TMAX,CGrt(oO)-CGM(65 J l/HTSS  +  1.01 
NREM  =  NHT 

29  IF  (NhEtM— NGIM  )  31,31,32 

21  NbTEP  =  NREM 
NkCM  *  0 

GO  TL  33 

32  NSTEP  =  NG1M 
KREM  =  NKfcM-NuiM 

33  00  3t  K*1,NSTEP 

TIME  *  HTS3*TLUATF iKHT)  +  CCM(b3J 
KiTfcP  =  K 
KH  I  *  KHT-*-  i 
CALL  GEGM 
CALL  GAVE 
3L  CALL  LUG  It 

^  -W  i 

ALBIS  OUTPUT 

•ML  =  1 
DUl  LT  *  1 
3 o2  LP  =  o 

Go  TL  (9oa ,4C2 »4c3 ,3U3 ,304 1 »KL 
**Ul  WRITE  uUTPoT  TAPE  b,8ul 
GO  TO  do  8 

<*1*2  WRITE  OUTPUT  TAPE  6,802 
Go  TC  303 

4u 2  WKlTL  OU1PCT  TAPE  0,803 
GU  Tb  30  3 

302  WRITE  UUTPul  TAPt  0,3yU/ SNAME, TNAME 
Gb  TC  d03 

WRITE  OUTPUT  TAPE  6,391,  SNAME,  TNArtE 

303  LL  *  u 
3d 8  KG  *  LT 

GU  Tb  (  8bl » 302  » 3c3 , 30b , 307  ) » ImL 

5bl  WRITE  OUTPUT  TAPE  o  ,9ul,  GT  (KG,  1)  ,  ( GT  (KG,  I  ) ,  l*?6»i  1 )  ,GT  ( KG ,  23 
1  GT (KG, 24  I 

GU  TC  do 6 

302  GPR( 1)  *  GT( KG, 2/ 

GPR(2)  =  GT(KG,4) 

GPR(3i  *  oT(KG,5) 

GPR(4)  *  GT(KG,3) 

GPR ( 3)  *  GT( KG, i*) 

GPKibl  «  GT(KG,13) 

GPw  ( 7 )  *  GT(KG,Lfc) 

Ob  35o2  1*1, 7 

33u2  GPK(i)  *  G PR (I)*37 .2937793 


r.  c 


,1  +  161, GUK 


bb  I  l,  : 
L  wK  1 1  l  u 
Gii  1L  J 

?  In  i\  1  f  h  U 
6  LI  =  LT 
IP  4  L  T-i\ 
S  LF  *  LP 
IF  (LF-<t 
l  LL  =  Ll 
IF  (LL-‘j 
i  W  K  a  T  t  G 
bb  IL  J 
c  KL  =  l\L 
IF(KL-i> 
H  C  C  i\  T  1 1\  L 
i  HJkMAI 


X  ^HK 
4HDE 


ikDQT  luX  5 
SfcC  7X  2HF 
9 X  6HFT/SE 


U  l\C  H  OVC'  tlO  ills?  ii  iOl  t 

ti^b  FUiNi-IaI  ( //(F6.  L  ,F9.1  ,PiU»i  )  ) 
0^1  L  U  IS  I  I N  U  t 

PLLlfUG  FuH  VtKSIUM  ll. 

iFlCCMM))  '»U,yu,7bu 
7b b  Ub  /Ik  K  =  l,iMSThF 

fllh/l  =  h(K,^)*1u., 


AZAAT 

IHDfcG  7X )  // 
ALPI2) 
3,2)  / 5HJ  S 
/) 

:  VB  S/N 
OELD  Dfc 


71  u  H(KfjZ)  =  (h(l\,0  2)+20.)/5. 

XHaX  =  CbM(to2)/oObu. 
utLTAA  =  XMAX/u«Ji 
^.cHliX  =  -UtLTAX 
XYAX(i)  =  U. 

X  Y  A  A  (2  )  =  C. 

YYAXllj  =  U. 

7  Y  A  X  ( 2  )  *  iG. 

XXAX(i)  =  XMAX 
XX AX  ( l  )  =  O. 

XPH ( i  )  =  -UcLTAX/lu. 

XPh(2)  =  oELTAX/10. 

OXi  -  -• l*UtLTAX 
DX2  =  -• 2*OtLTAX 
DX2  =  -.a*l)ELTAX 
UX*<  =  -.4*ocLTAX 
DX  5  =  -•  3*DELTAX 
MPY  =  uuM(o2)/C0M(Od  +  .Oi 

CALL  PL  I N1  (-l,6HUWHlTt,-6,3,0X4, 1 ,  Z  EKCJX  f  DEL  TAX ,  XMAX ,  9.9,  l,3Hl.O) 
CALL  PLALP  (-0, 2,0X3, 7. 9, 2H.b) 

CALL  PLALP  (-0,2,uXJ,:>.9,2H.6) 

CALL  PLALP  <-o, 2, 0X3, 3. 9, 2H. 4) 

CALL  PLALP  (-0,2tOX3,1.9,2H.2) 

CALL  PLALP  (-6, 1,0X1, -di  1H0) 

CALL  PLFUN  ( 2,U,XYAX,2,YYAX,1, 1) 

Uu  72c  1  =  1 , 1C 

YPF(i)  =  lO.-FLUATF(I-l) 

YPH  (2 )  =  YPH(l) 

72C  CALL  PlPUN  ( 2, 0 , XPH, 2 , YPH, 1 , 1 ) 

00  7c3  1  =  1, NFV 

XPV(l)  =  CCM(6l)*FLUATF( 1 ) /6ubu. 

XP V ( 2 )  =  XPV(l) 
lF(XMUDFd,5)  )  721,722,721 
7*1  YPV(i)  =  -.03 
GO  TC  723 

7 22  TEMP  =  FPbCDF(XPV) 

YPV(l)  =  -.1 

CALL  PLAlP  (-4,o,XPV+UX3,-.23, TtMP) 

723  YPV(2)  =  -YPV( 1) 

CALL  PLFON  1  2,0,XPV,2,YPV,  Id) 

723  CUNT IiMCE 

YXAX(l)  =  0. 

Y  XAX(2 )  =  0. 

CALL  PLFUN ( 2 ,0 , XXAX, 2 , YX AX , 1 , 1 ) 

CALL  PlFUN  (2,0»Hd,l)»NST£P,H(l,2),l,l? 

CALL  PLtNO  (1) 

CALL  PL1N1  1-1 , CHOrtHl TE,— 6 , 2 ,0X3 , 1 »ZEROX, OELT AX, XMAX ,9.9, 1, 2 H 30) 
CALL  PLALP  (— b,2,0X3,7.9,2H20) 

CALL  PLALP  ( -0,2 ,0X3 , 5 . 9 , 2H10 ) 

CALL  PLALP  (-0,1, 0X2, 3. 9, 1H0) 

CALL  PLALP  (—0,3, 0X4, 1.9, 3H-10) 

CALL  PLALP  ( -0,3, 0X4,- .1 ,3H— 20) 

CALL  PLFUN  ( 2,0, XYAX, 2 ,YYAX, 1, 1) 

DO  73C  1=1,11 

YPK1)  =  1C.-FLUATF(  1-1) 

YPh( 2 )  «  YPH(l) 

730  CALL  PLFUN  ( 2, 0, XPH, 2 , YPH, 1 , 1) 

DO  735  1=1, NPV 

XPV(l)  =  COM (Ol ) *FLOATF l I ) /oubO. 


AP  V  U  )  =  XPVU) 

IF (XFUUF (1*5))  731.732.731 
7  31  YPVU)  =  3.95 
Gd  Tl,  735 

73 ^  TcFP  =  FPbCOF(XPV) 

YPVU)  =  5.9 

CALL  PLALP  (-9,o,XPV  +  UX:>»3.75»TtMP) 

7 j3  YPVU)  =  b.-YPV(l) 

-ALL  PLFulNI  U,U»XPV»2»YPV»1»1) 

735  CCM1NUE 

YXAXU)  =  4. 

YXAXU)  =  4. 

CALL  PLFUN  U,U,XAAX,2  »YXAX,  1,  1) 

CALL  PLFUN  ( 2* C , Ht 1, 1) ,NST £ P ,H 4 1 ,32 ) . 1 , 1 ) 
CALL  PLENU  IKtND) 


9C  IF  ( hi  K  L  i"1 )  9,9,91 
91  IF (CCrtlOi)  )  29,29,19 


l.i 

SU/bKUU  1 1  (\|4_  TAROIM  TMAX) 

DIMENSION  TNAMEI2) 

DIMENSION  EPHI 322,6)  yCROStCI l HOI » 2 ) » GT ( BO , 32 ) , H( 80 , 40 ) , COM( 1G0) 
COMMON  TIME,  K3TEP,KHT,FPH, CRUS  EC,  GT,H,COM 
C  TARGET  INPOT 
luUl  HUP MAT  ( 6X , 2 Ao , 2X, I3y2XyiF3.0) 

10U2  FORMAT  < I3,3X,1F4.G,2<3X,1F7.G,4X,1F7.1,2X,1F6.2) ) 

2ot  1  FORMAT  ( 1H1  ,2Ao,24H  bUOST  TRAJECTORY  COORDINATES,//) 

20u2  FORMAT  ( IH  ,24X,10HDUWN  RANG Ey 2oX , 8HALT I TUDfc , // ) 

2oui  FURMAT  (iH  , 10HCARD  T IME , 2 (4X ,3 1HPUSI T ION  SPEED  ACCELERATION, 

1),//) 

2uo4  FURMAT  1 1H  y  14, 2X,  F4.0 ,2  <4X  ,F8 . 1 , 2X,  F7.  i ,  3X,  F6. 2 , 4X)  ) 

C 

TMAX  =  w. 

10  Kb  AL.  InPUJ  TAPE  5,1001, (TNAMfc(J) ,J=1, 2) ,NMAX,TAU 
WRlTc  uUTPoT  TAPE  6,2uol ,  (  TMMEI  J )  , J«1 ,2 ) 

WRITE  uUTPUT  TAPE  6,2u U2 

WRITE  OOTPOT  TAPE  o,2uu3 
DU  11  K= 1 , NMAX 

READ  IimPUT  TAPE  5 ,  iuC>2  ,L ,  T ,  I  EPH(  L ,  J ) ,  J»1 , 6 ) 

TMAX  =  MAX1F (T ,TMAX) 

11  aRITE  OUTPUT  1 APt  b, 2004, L»T,(EPH(LyJ),J= 1,6) 

CUM  I  4 )  =  TAJ 

RETURN 

ENC 


1.2 

SUBROUTINE  CKOS  I N 

DIMENSION  EPh(  322tO)  ,tKUSEC(lt!01,2)yGT<80,32)  ,  HI  80 ,40  ),  COM!  100  ) 
COPMLN  TIMc, KST tP ,KHT , EPH, CHCSEC ,  GT , M, CUM 
C  CKUSO-OECI 1U  INPUT 

loOS  FURMAT  (13F0.i) 
lout  FORMAT  ( IF 0 • 3 , 12Ao ) 

4U0C  FORMAT  I  IH  yU<F6.1,AA),2X) 

4Uo2  Format  UHl,l2Ao) 

OlMcNSiUN  EX ( 1 3 ) yKNAME ( 12) 

C 

30  READ  INPUT  TAPc  3, ioOo ,DELALP, (KNAMt (Ki , K* 1 , 12 ) 

WKilE  OUTPUT  TAPE  o,<KW ,  <  KNAMEC  L ) , L*1 » 12) 

K-U 

31  READ  INPUT  tape  3,iuu3,tEX<J),J=i,U) 

Ir(cX(li-loa.)  32,34,34 

22  WRITE  OUTPUT  TAPE  6,4uoo,(EX(L)yL=l,13) 

UU  3  3  J=  1 y  6 
LH  =  4.»J 
LV=LF+1 
Lt*K+J 

CKCStU ( Lb , i )  =  EXPFI IcXI Lh ) /10. ) *2 . 3025631 )*iu . 76387 

33  CR0SLC(LL,2)  =  EXPFI  ( EXI  LV )  /10.  )  *2 .3023631 ) *10 . 76387 
K  =  K  frc 

GO  TC  jL 

34  UON  J 1NUE 

CKCSEu I K , 1 )  =  CKCSEu  ( IS”!  ,  1  ) 

CRC SLU I K , 2 )  *  ChuStC(K-l,2) 

UOP(ii)  =  OcLALP 

Kh 1  URN 

END 


-  90 


o  r>  n 


SUcnLollNb  UboM 

UlMof.SloN  KobI{  3)  ,  VuoT  (3)  ,KttM3)  ,  Vu  Ja  (  j  )  ,  uRGc  H  3)  ,UV0JT(  3) 

01 Kl No  ion  rtL  AT  (  j>)  ,  RUT  A(3  )  ,  UkUA  J  (  J  )  , u VUaT  (  3  )  ,  UKOTA  (  3) , UVUUA (  3  I 
ulFcNSlon  VLUbK(i)  ,UVtKTt3)  ,  I  K I  AuA(  3, 3)  ,  HAT  (  3  )  ,  hAB(  3  )  ,R  00b  (  3) 

U  1  KfcnSiuN  HULL  (  3)  , KuAo ( 3  )  ,  RLmC  (3  )  ,  kGTo (  3)  •  UK  G  All  (  3)  ,UKUAC ( J ) 

Dl  P  ti\3  iuN  LKoI b l j )  »ALP(j>),TFblA(;>),  S  lOCkCi  (3(2) 
ui Me Nil  UN  VoaT  (  j  )  f  kuTC  (3  )  f  Uk(JTC(  3  ) 

UlPcNiluN  VbCTl( 3) ,VboT2< 3) 
blPfcNSiUN  HALO) 

blPtRSiuN  Lr-nt  A<.f  a)  ,CkUSbC(  ibL  i*  4)  ,  oT  ( tio  (  32  )  »HC0M(ac  (40)  » C  UM(  100 ) 
LbMMLN  T  IRE »  kS  I  eP  »  KrlT  ,  cPH »  LKCSLC  ,  bT  »  HCOM ,  COM 
C 

C  Co  A  S  1  AN  I  S 

C 

K  =  KSlbP 
UVbkT (  l)=o. 

UVbk  T  ( 2)  =C. 

UVbKT(  3)  =1. 
ke  -£.76  7b4oC  • 

PI =4.*ATANF( I. ) 
kAlLCb=16o./Pi 
CLloHT=u.Vo357ioOb+u4 
VV  A  =  CoM ( 2 ) 

OLLALP  =  CCM(ll) 

C 

OtCMb  Tk Y 

TAKOtT  ANb  AlkCkAFT 

CALL  TAkOcT  ( T IME, ROOT , VOGT) 

CALL  AIKLFT  (  T  i  fit  t  RuOA  » VGUA ) 

CALL  UNIT  (ROUT, GROUT, RUT) 

CALL  UNIT  ( VLUT,UVbUT,VVT) 

IF  (  T  i i“l t3 )  2  24,22  3(224 
223  DU  223  J  =  l,3 

UVCOT(J)  =  UVtRT(J) 

223  UKCuT(J)  =  UVLKT(J) 
ciH  CONTINUE 

00  leO  J  =  l,3 
KuAT ( J ) =RUGT ( J )-KUGA( J ) 

RGTA( J )=-RCAT ( J  ) 

22C  VUAI(J)=VOLT (J)-VUOA(J) 

CALL  UNIT  <RUaT,URCAT,RAT) 

CALL  UNIT  (VUAT,UVGAT,VAT) 

CALL  UNIT  (RCTA, UKOTA, SCALAR) 

CALL  UNIT  (V0uA»UVG0A»VVA) 

CALL  CROSS  < LVEkT, UVUUA, VECTOR) 

CALL  UNIT  (VELTUk, VECTOR, SCALAR) 

DO  221  J=i,3 
TkIALA(1,J)=UVUUA(J) 

221  TRIALA(2»J)=vECTCR(o) 

CALL  CROSS  ( 'JVUUA,  VtCTUR, VECTOR) 

CALL  UNIT  (VEC TOR, VECTOR, SCALAR) 

DO  lei  J=l,3 

til  TKlADA(3,J)=VcCTCR(J) 

CALL  TRIAD  ( TRIADA,UKUAT,ELAAT,AZAAT) 


c 

c 


BUUNCu  P  C 1 N  T  -  kLUnD  CAKTH  WITH  REFRACT  I  ON 


<S=3wkTKRuAT  (  i)**<i+ROAT(2  )**2i 
A=  l  • 

B=-1.3*kS 

C=-KE*(kuCl ( 3)+KLuA( 3) )+u. 5*RS*RS 
D*Rfc*KS*KUCA(:>.» 

Kb  1  =  (RUUA( 3J*k$)/(RuuT (3 i  +  RLUA (3) ) 
cjI  F=kbl*(RSi*(RSl*A+o)+C]+U 
PP=RSi#(KSi*3.*A+*  .#d)  +C 
UELRS1=-F/FP 

IP  UfcSF(ULLKSi  )-i.  )  tii *232, 23 < 
c32  K3i  =  RSl  +  L)ELRSi 
GO  i L  2j1 

c3j  ALPHAd*KOUA(  3)/RSl-KSi/(2.*RE) 

H*RS*kb/{2.*iu:] 

Hl=KSl*KSi/(2.*kt) 
kS2  =  k3**RSl 

K=  Swk  IX  i  kuGmi  3 I-kuJT  { 3)  +h  )**2+ks**2 ) 

Ki=bwkTP( ( KUCA(3)+Hi ) **2+kSi**2 ) 

R2=SGRTP ( (ROOT ( J**2*RS2**2) 

ALPHAC  =  ARCSIN( KUJA(3)/RAT  -  ( RAT**2-RuUA< 3 ) **2 ) / ( 2. *RAT*kE H 
C 

C  TIKE  ucK  I VAI IVES  uP  R3»K  »R1  AND  R2 
C 

K SKSuT  =K (JaT (  I  *  #VGA T  ( 1 )  +RUAT  ( 2  I  *Vl)AT  ( 2  » 

kSCCI =KSKbOT /KS 

BuCT=-i.5*KSUuT 

C UU  =— Rfc  * (  VULiT  ( 3 )  P VUOA  ( j  H  «-RSRSDT 
DDLi  «ist*(K3ULr*RUjA(  31+RS* VOGA (3 )  ) 

KSiDUl=-(KSl*(KM*bDUT+CDUT)«-QCUT>/(RSl*(R$i*3.*A+2.*B)+C) 

•  k32uUT=kSUUT-KSlbUT 

RKUuT=(RuGA(  3)-KuUn  3)+H)*(VUUA(  3 } -V0UT{  3  )  «•(  RSkSDT/R  E  )  i  ♦J’.SRSDT 
KuCT*kRDUT/K 

KikiCT-(  RUGA  (3)+HjL)*(VGGA(3)  +  (  RS i*RSiOUT  /RE )  ) +RS1*RS1D0T 
RiCU  J  =  R1kJlDI  /KI 

K2K2D  T*(  ivbaT  (  3  J-H+.UJ  *  ( VOuT  ( 3 1  -K  SRSDT/RE  +  *  RS1*RS1D0T/RE )  )  + 

1  kS2+RS2DUT 

R2CGT=is2R2DT/R2 
TOlK=(2.*k)/CLlGhT 
TUbUR-  (2  .*(R1+R2I J/CLIGHT 
T  S  ING=v>«  3*  { TCIK+TOGUb) 

OELT=  TGUUb-T  SING 
UUP01R  =  2  .*ROCT 
UuPi^rtUUT+RUUT+R^DOT 
DOP2=<..*(RlDbT+R^DUT) 

C 

C  BuLNCE  AND  CLUTTER  PUINTS  IN  PLAT  EARTH  COORDINATES 
C 

HaT(  1)=KUgT( l)-ROOA(i) 

HA  T ( 2) =RUOT ( 2) -KCUA( 2 ) 

HAT(2)=U. 

CALL  UNIT  (HAT  fHAT  f SCALAR! 

CALL  MULT(RSltHAlfHAb) 

ROCfcUi)=RGGA(l)+HAB(l) 
kUUU  ( 2  )=K(JUA (  2 )+HAB(  2 ) 

RUcd(j)=U. 

Xc  *  SURTF ( RA  T**2  -  RGGA(3)**2) 

CALL  NUL  T ( XL  »hAT  #HAC I 


RulC  (  i  )  =  RULA  ( 1  t  *  HAc(  1) 

RULcU)  =  KCGrt(Z)  ♦  HAL ( 2 ) 
kuLc ( j ) =0 . 

DU  2Ac  J  =  1  » 3 

kUABI J) =KLuB ( J )-RUuA( J ) 

ROACI J )=KUUC ( J l-RUUAI J) 
kJ  IC(J)  =  kLUH  Ji-KUUTl  J) 

2ac  Kulbt J )=kLGb( J )-RUUT( J l 

UALc  UnIT  (KuAb»  Ut\UAB»KAb  ) 

Call  unit  (alac» lrgal »rac) 

Call  unIT  (RUTb,UKCTB»kTB) 
call  CHIT  <RUTC,UkUTC*RTC) 

CALL  fklAU  (  TklAUA,URUAB,ELAAB,AZAAb) 

CALL  TRIAL  (  Tk  I  aUA  «  UnU  AC  « E  LaAC  » AZAAC  ) 

C 

C  DIRlcT,  SlNGLc.  AND  UuUbLE  BOUNCE  SIGNALS 
C 

Call  oci  ( uvcu; »uruta»CuSalp) 

CALL  CnUSS  ( UVCCT , UKUT A, VtC TUK ) 

CALL  UNIT  (VECTUR»VECIUK»SINALP) 

CALL  ANGLE  { SlNAlP,CUSALP,ALP( i) ) 

LALL  CUT  (UVCUl  »  UKU  TtJ fCuSALP  ) 

Call  CkuSS  ( UVUcT»UKUTo»VECTCk ) 

CALL  UNIT  ( V LCTUKf VEC TUR » S INALP ) 

CmLL  ANGLE  I SiNALP  »CUS ALPi ALP ( 3} ) 

1MTIME)  2A2,e41,2a2 
cAi  ALP  I  3 )  =  AL‘Mi) 

2a  2  Cuk  I  INUl 

CALL  CkUSS  ( LVuuT » KOTA, VECT1 ) 

CALL  CRuSS  (  cV EK T » RUT A » V EC T2 ) 

CALL  UNIT  (VluTI  ,VECTi , SCALAR! 

CALL  UNIT  ( VeCT  2  » VECT2 , SCAl AK I 
CALL  CUT  ( VtC T  I » vcCT2 »COSTh ) 

CALL  C.<U  Si>  I  VECfl,  VtCT2,  VEC1CR) 

CALL  UNIT  ( V EC Tuk » VtCTUk t S INTH ) 

CALL  ANGLE  ( SINTHi cUSTHi THeT j 
CU  2a. ^  J  =  1  »3 
c  A  S  ThcIA(J)=lHET 

ALP  l  c  )=ALPl 1 )  +  U. 3* ( ALP ( b )- ALP (  i)  ) 

Uu  2  bu  J  =  1 1 3 
ALPP=mLP(J)*RATuDE 

CALl  inaTSuT  (ALPPfSlGMAHf  SIoMAVtUELALP) 
a  1 NTH  =  s I NF  (  THeTAIJ  )  ) 

CuSTE=cuSl_(  TncTAlJJJ 

SI  CcRut  J  «  i.  i  =  SiGMAri^cUSTrl^^  +  S  1  GHA V*S INTH**2 
Lbc  SiLCKUlCf  2)  =  SIoNAH*SlniH**Z  +  SIGi-1AV*CGSTH**2 
C 

C  STlacC  uecMEIkIc  VARiAbLES 

C 

GT(Kfil=1iME 
GT  (  R  t-i  )  =  tLAAT 
GI(K,2)  =  AZAAT 
GT (K»a  5  =  eLaAo 
GT{K»j)  =  tLAAC 
GT  ( is  »u  )  -  R 
Gflis,/!  =  R I 
ol (K  »b  I  -  n2 
GllKt H  )  =  kCL  T 
GT(Kfic)  —  k  1 1 Cl 
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oT  (KfilJ  =  K2LIGT 
L>T  tis»  1^:J  *  ALPhAB 
G1  IK, 13)  =  ALPHAL 
bllKn-t)  =  ALP  ( 1  i 
Gl  (K,15)  =  ALP ( 2 ) 
GTtk,io)  =  ALP(^J 
G1(K,*7)  -  SlGLKLil,l) 
GT(k,1cj)  =  S1oCku(2,1) 
bTCktlv)  =  6IoLi\C<3,i) 
GT (K,iol  =  SloCkLt 1,2) 
GT  IK ,£ 1 )  =  SiGCRC(2,2) 
GT(k,c2)  =  SIGLkC( 3,2) 
GT (K  ,23 )  =  TCik 
GT  <K,2*«)  *  ULLT 
GUk,25)  *  ULPOI  i< 

GT Ik  ,2c)  =  DCP1 
Gl  (K,27)  =  0LP2 
gi  ( k , 4 o j  =  rntTAtu 
GT(K,2S)  =  KLLT ( 3 ) 

Rt lukw 

cNC 


SUBROUTINE  TAkotT  (TIMt,R,V) 

UlMtkSlON  H ( 3) ,V (3) 

UlPENSiUN  tPhl 322,6) ,CRGSfcG< ibOi, 2) , GT( 30 , 32 ) ,H( do , 40 ) ,CGM( 100 ) 
CUPMLKi  UUMT,K3TEP,KHT,tPM,CRUS£CtGT,H,CCJM 
TAL  =  GUM(<J) 

K  *=  1  iMt/TAU  *  1.01 
Tl=TAU*PLGATf <K-i) 

TG*1  IMt-II 

Kill  *  tPhlK ,1 ) ♦LGM( 35) 

M2)  =  LPH ( K , i ) *CuM( 34 ) 
k(2)  =  cPIKK  »<t ) 

.  V(i)  *  lPH(K,2)<iCGM(  ) 

V  (  2  )  *  t-PH(K,2)*Cu(*(34) 

V ( 2 )  *  LPhlK.a) 

1 f (CLM( 6 1 )  )  3,3,2 
3  INTO  <.,2,1 

1  ThL  IA=  fG/T  All 
Tht  TA2= I HlTa**2 
THfcTA3*ThLTA2*THtTA 
A=2'.*THtlA3-3.*THtTA2'»'l. 
b=-2.*Thfcl  A3+3.*TH«iTA2 
c=TaG*<  rhtlA2-2.*THETA2+THcTA) 

U=TAU<t(THtTA3-ThtTA2) 

Ki=K+i 

M2)  =  A*fcPH(K,4)+to*EPH<Ki,4)*C*tPH(K,i>)+0*tPH(Kl,5) 

V(i)  =  A*‘tPH(K,S>')+8*fcPk(M,5)  +  C*EPH(K,6)+0*tPH(Kl,6) 

Xk  =  A«‘LPH(K,i)+o*tPH(Kl,l)+C*fcPH(Ki2)+l)*tPH(Kl,2) 

XV  =  A*fcPH<K,2)+b*EPH(Kl,2)+C*fcPH(K,3)+U*fcPMKl,3) 

M  1)  =  XR*LGM3i>) 

M2)  =  Xk*LUM<34) 

VC1)  =  XV*CoM35) 

VU)  =  XV*LUP(j4) 

2  RETURN 
tNL 
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OUcisLU  T  i  Nt  AlKCbT  (flME,K,V) 

L  1  H  b  b  b  1  O  i\  k  (  2  I  *  V  1 3 ) 

D  i  k  LNb I uN  bPh(3^2',6),CKJSLC(lbul,2),GT(8u,32)  *H(du,Au)  ,  COM(  100  I 
CuRMLN  DUMT,  KjTtP,KriT  ,  bPH  ,  ChCbbC  ,  0T ,  ri ,  CUM 
V  (  A  i  =  CUMjl) 

Vic)  =  CuM(3Z) 

V(i»  =  CUM33) 

ThLTA  =  CLM(43)*( TINE-CUMI bi » 

Kli)  =  COM(oa)  +  lritTA«V(l) 

K ( <: )  =  COK(t>)  + JhtTA*V(  Z) 

K ( 3 )  =  ouM(7H-lhbTA*V(;>) 

Kb  7 Uki\ 

ENL 


SubKLC  1 1  i\t  UNIT  (XX,YYVZZZ) 

Gif'  bNb  IuN  XX  (  3  )  , YY (3 ) 
ZZ2=SCNTF(XX(l)**2+XX{2J**2+XX(3f **2) 
IFlcli.)  1,1,2 
i  YYli)  =  U. 

YY ( 2 )  =  0. 

Y Y  ( 3  )  =  U. 

(30  TO  3 

<.  YYli)  =  XX( 1  \/ LLL 
YY  (2  )  =  XX ( 2  i/ ILL 
YY  (3 ) =XX( 3  )/ lLL 
3  RtlUKN 
tND 


bULRGCTlNt  C*USS  (XX,YY,ZZ) 

oi Pension  xxij) ,yy(3) ,zz(3) 

A  =  XXt2)*YYlj)  -  YY (2J*XX( 3) 
b  =  XX(3)*YY(1.»  -  YY(3)*XX(1) 

ZZ  ( 3  )  =  XXm*YY<2)  -  YYU)*XX(2I 

zzm  =  A 

ZZ (2 )  =  B 

kb  TORN 

bNC 


SuuKUuT I Nb  TRIAD  < TK l , VEC , EL » AZ) 
DIMENSION  TRi (3,3) ,VbC(3) ,UVEC43) 
CALL  UNIT  ( VbC ,UVbC, SCALAR ) 

SuM  =  U  . 

SUR2=u. 

SGKj=U. 

OU  1  J  =  1 » 3 

SURi=SUMi+TR141,J)*UVEC(JI 
3JM2  =  3oM2  +  TRI(2,J)*UVEC(  J) 

1  SUM3=SuM3+TR 1  (3,J)*UVbCiJ) 
CUSLL=bGRTHSUMi**2+SUM2**2) 

S I NAZ-SUH2/CCSEL 
CU  SAZ  =  SUM1/CC$EL 
bL  =  A  Tank  S0iM3/CoS£L) 

CALL  ANGLE  C SINAZ.COSAZ, AZ I 

RE1URN 

ENC 
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o  rr> 


SUEkUoTINt  MCcT  (XXX,YY,ZZ) 
OlPENblUN  YY(.>),ZZ(3) 

ZZ ( 1 )=XXX*YY ( i ) 

ZZ(c)  =  XAX*YY(<.) 

ZZ (3 )=XXA*YY (j) 

Rt 1UKK 


SubRCUl  lNt  OCT  ( XX  » YY  » ZZZ  ) 

UlPtNi»iuN  XX 1 3) » YY (3  ) 

zZZ  =  XXU  )*YY{  i)+XX<2)*YY(2)+XX(3)*YY(3) 

h£ IURN 


SubKbU  I  I  Nt  AMiLt  IS1NCHI  ,C0$CM ,CHl ) 

PI  *  3.1415927 

R=5URT  f  (  SiNCHI**2+CG$CHl  **2  ) 

S 1  NCR  I  =  S INCH1/R 
CuSth 1 =CCSCHI/K 
IP  (CCSCHi  )  S,l,4 

1  lP(5l*CHl)  2.2*3 

2  CHI  -  (  3 •  * P I  J/2. 

Gu  (L  10 

2  Chi  =  PI/2. 

Gli  1C  IC 

4  CHIP«A  TANp  ( SlNCHl/COSCrll ) 

UMSlNChl)  8,5,5 

5  IMcoSCHl)  7,0, o 
t  CH  1  =  CH IP 

go  rc  iu 

7  CHI«CH1P-»PI 
Gu  rc  lo 
iPiCCSCHl)  7,9,9 
CHI  *  CHIP+Z.+Pl 
it  CONTINUE 
RETURN 


GUbKUUTINL  RATSCT  (ANGLE  »SIGPAH,  S1GMAV. OELALP) 

OIMt.NSJ.UN  tPh(  322f  6)  »CKOSEC(  I8C 1 ,  2) » GT(  80  ,32 )  ,H(80,40)  ,COM(  100) 
CUMMcN  T IHt»KSTtP.KHT,EPH,CRCSfcC»GT,H,COM 
K 1  =  ANoLt/UtlALP  ♦  1,5 
SIGMAh  =  CRGSEC (Kill) 
iIGMAV  =  CKOStC(Ki,2) 

RETURN 


v  ■hi -vmp 


.,  (waW  *t*WWrtO»*W*il»4tW,»^«v^MtiWiSS',<wij,  *,tfj,rs»  i:mm»/t  MeWMW  o-uiWI.E. 


1 


H 

0 

0 

n 

o 

o 

o 

Q 

D 

Q 

a 

o 

3 

3 

| 


SUUKUOTINC  U AVC 
UlMCNSlON  TPK (  12  ) 

uiMthSlON  tP H  (  J^to)  y  C  k(j  SOC  ( 1601 y  2 ) y  GT(  do  iic)  yH(80  y  40  )  y  COM  ( IOC  ) 
COMMON  TIMtyKSTtPyKHTyhPHyCRUSECyGT  y  HyCOM 
K  =  K5 1  CP 


C 

C 

C 


c 

c 

c 


GtT  inputs  from  common 


VVA  =  oUM(2) 

FMC  *  CoM  ( 12  ) 

HAPtK  =  CGM(  13) 

VaP  tK  =  COM ( 14 ) 

AAPtFF  =  CGM ( 15 ) 
tTA  =  COM  (16) 

PCl\  =  COM ( 2u  ) 

PPtAK  =  COM ( 23) 

FK  =  COMU4) 

T  u  Ik  e  L  L  =  COM (25) 

SYSLF  =  CoM{ 2o ) 

TZbtAM  =  COM  ( 27  ) 

FNL1SC  =  ouM(2b) 

HT  Si  *  CUM<2S) 

VSC  =  COM (30) 

XlHUtL  =  CCM (37 ) 

XiVUfcL  =  COM ( 3d ) 

XIFZ  =  C0M(3y) 

X1V2  =  CoMUGJ 
XIF2  *  oOM (41) 

XiW  *  CUM (42  ) 

StA  =  COM ( A3  ) 

POL  =  CUM(44) 

FNTRK  =  CUM(AO) 

CRL  =  CUM(47) 

bN Sk  =  SsikTF  (C0M(4G)  ) 

FnS  =  OuM(A<>) 
bi  =  COM (bo) 

ATlIS  =  CCM(3x)**2 
AlEPS  =  CoM( 52 ) 

ATtoS  =  COM (33 ) **2 
S i GOP  1  =  COM ( 5o ) 

S1G2  =  10.*>MCCM37)/iU.  i 
RliLANK  =  COM  ( 64  ) 

RAiN  =  CUM(2i) 

SUMS  iG  =  lU«<tJMcCM(22)/10») 

GtT  STUFF  FROM  GcOM 


L  =  K 

KANGt  =  GT (L y6 ) 
Kl  =  G  i  (  L  y  7 ) 

R2  =  GT(L,d) 
IF(POL)  14yi4,lA 
12  S1C0C  =  G1 ( L  y 1 7 ) 
S  i  Gob  =  G1 (Lyidl 
SiGbb  =  GT(LyiS) 
GO  TO  10 


1 A  S  l GuC  = 


-  f' 
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i 


GT (L  y2o) 


rr  r 


MkUHlWi 


S 1  G  J  e  s  Gl(L»Ail 
blued  =  oT(L»22) 
ic  ALFhAc  =  bTlLtiJ) 

ALFmau  =  Abif  (G  I  (L  ,  12») 
XI  vc  =  gT(l,:jJ 
ouUo it  -  gtil.zb) 

uUFi  =  G  T l L , 2o ) 

OOP,:  =  bT(Lu7) 
tLAAT  =  Gi (Ltcl 
tLAAU  =  G T ( L  t  A ) 
uLlT  =  GT(L,2A) 

C  =  S.ojb71  td 
PI  =  3.1415*2/ 

bGLlZ  =  i.dduA  c  — 23 
SUMP  =  290. 


ANltt\NA 


C 

c 


PLA.-lti  *  o/(FMC*l.t6) 

ANAKEa  =  AAPtFF*,HAPc.K#VAPfcK 
AuAlUA  =  A.*Pl*Ai\AHcA/(FLAMti**2) 

Ht  (Ah  =  TLAMb/HAPfcK 
bCTAV  =  FLAMb/  VAPeK 
HbfcAMP  =  1. 

VbLAMP  =  1. 

IP  UBS*  ULAAT-tlAA6)-bi£TAV/2.)  5ul,5dl,:>u2 
5wl  liucAMP  =  'i. 

Gu  TC  Bui 
B02  odfcAMP  =  tTA 
5u3  Cui\l  livUt 


FSCli  =  AlhOtL**!  VOcL/<FlUS*UtTAVi 
I<-i*P  =  PbCHMTlrtfc-TZbcAM)  ♦  XlHZ 
XIFG  =  TfcMP  -  XI HJtL*l NT FlTtMP /XI H2) 
1F1VSCJ  32u»5lci 52u 
510  XlvG  =  XIVZ 
GU  TO  3iij 


Bee 


BBC 


ftFP  =  XIVZ  +  btTAV*INTF(<FSCh*( TIME- 
XiVG  =  TtMP  -  XlVUtL*INTF( TbMP/XI V2) 
CCNTlNUc 

FUP  NGw,  life  X1HT,  XIVT 
XIFG  =  GT(ltB) 


TZB£AM)+XIHZ)/XIH2) 


XiVG  =  G  f  (  L, 2  ) 


Phil  VC  =  XlVC-XiVG 

c 

C  KtFLfccnuN  CULFF iCItNT 

c 

1F1PLU  72c,  7jLCf  720 
71c  IF  (ALPhAb-i.GA (2)  711,711,712 

711  KtFl  =  1  .-ACPHAB*,6/P1 
GU  TL  / 30 

712  ReP'L  =  .8 
GO  Tu  730 

72C  IF(ALPhAb-.Ob727)  7^1,721,722 

721  KfcFt.  =  l.-ALPHAb*32.4/PI 
GO  TC  730 

722  kLFL  *  . 8-.V2bi6*tXPF  (-10.  O8*ALPFIA0/PI ) 
73t  IF (SEA)  73t, 731,732 


p>  ■ 

'*  •? .  i\ 
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ri  n  o  r>  r  n  r  nr  <r  r>  r- 


73  i  TtPP  =  i. 

Go)  T  L  7  j>  J 
722  TtWP  =  1  ti « 

/3J  ntFL  =  K4Pl*LXPF(-o.#<Pl*TfcMP*ALPHAB/PLAMB)*#2) 

MOnAL  PKLCfciilNG 
Pci  =  PON 

MlPcSPHtRlC  ATTfcNclAT luN 

ATMOU  =  1. 

ATMC3  =  i. 

ATFlcu  =  1. 

TkAwiMlTTLK 

PZ  =  (PPtAK*ANAKtA*#2i/(4.*PI*FlArtB**2*HANGt**4) 

Pa  Vo  =  PPfcA*.*TOktU.#FR 
PuLbt  *  T  G  Wfc  LL  /  PCA 

ilCisAL 

FLLii  =  3YSlF*P0S 
TfcPP  =  PZ*FcC:>S 

SQL  -  RtoP#S  lG0b*ATMlUU*HbfcAMP***.*VBEAMP**2 
i>UB  =  RMP*SlGUb*(RArtufc/(ki+K2) )**2*RfcFt**2*ATMLDB*HBEAMP**2 

1  ♦ViitAMP^BBLAMP 

ibb  =  ltMP*i>lGGB*<KANGE/(klFR2))**4#RfcFl.,J,*4*ATiMLB8*BBfcAMP**2 
IF (ucL  T-PuLSfc )  21J,2io,22o 
410  IF(ihA)  440,2 11,220 
211  ST1  *  5 jO+S6B*4.*SOo 
Gbl  ( C  430 
240  3IT  =  SOU 
2 3 u  CliMlNUfc 
C  » 

c  oliPPLbR 

C 

fUL  -  COPOlR/FLAMii 
Fub  =  UUPi/FLAMb 
~  -  Fbfc  =  GGP2/FLAMB 
0 

C  PFiASfc  c  INFERENCE 

c 

PHA  *  (K1+K2-RAi\0c1/F0AMB 
IF(PLL)  451, 432, 251 
2si  IF (ALPFAB-.ob727)  232,252,253 
452  PhA  *  PHA+.3 

chi  PlIAuEL  =  (0UP2-uCPi)*bETArt/(FtAMB*FSCHi 
C 

c  UUlt« 

C 

1 F  (  aLPFiAO  )  254,254,295 

254  SIGZ  =  u. 

GO  TL  332 

255  IF (AlPHAc-PI/2, J  247,257,29o 
25fc  ALPhAC  =  Pi-ALPHAG 

24  I  1 F (AbiF ( ALPHAC-P 1/2. )-l. t-5 )  296,299,299 
44B  AL  PH AO  =  PI/2.-1.E-3 
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tjM 

■ 


C 

c 

c 


c 

c 

c 


c 

c 

c 


2-iS 
-Sc  1 

3u2 

303 

Jo4 

305 

Jut 

Jut 

310 


Jc 1 »  ,?G2 , 333 


Jit 


330 
33  i 


COM  iNut 
IHSIouPT ) 

0 1 . 0  f  T  K  =  u. 

Co  To  330 

1  F  < IF Mo-3oo 0»  J  joc  ,  jui ,  _>y  J 

1  i~  (  S  tr  A  J  3u5, 304, 

bStA  =  .3 

OU  TO  Jco 

tJitA  =  3„ 

SIC^  =  2.*SlKF<ALPHAU*l0.**(~3.2  +  .b*aSEAI 

w  It  Sj<: 

lHSfcAJ  3cC,31o,320 
01  =  • woSo  1^ 

C<.  =  2.98 

=  •  CO 

It  3 JO 
=  •  lt>bo 
*  2.ob 
4.27 

*  C1*<AU>HAC**CM/<FLAMB*«<C3> 

=  KAivGc*C*rOwtLL/(^.*COSF(ALPHAC)l 

TcMP*bETAH 
TtMP»PI 

*  S1GZ*mCM 

*  S1GZ*ACS 
PZ*  S1GCM*;>YSLF 
P  Z*S1GCS*ETa**2*SYSLF 

CIS 


oJ 
GO 

Ci  = 

c*  - 

C  3  = 

Si  GZ 
TEMP 
AG* 
ACS  = 
SlGoR 

Sloes 

Cl*  = 
CIS  = 

cut  Tit 


->38 

340 

330 


C 

c 

c 


iPlRdlANK)  J50,33d,35u 

LLL^f  <(’,M,‘,'!S!'BeT4V/2>)  34O.3S0.3S0 
rLh!TK  *  CLUTTh  +  ClM  -  £TA**2*CIM 

*  Cfc  fAri*SIONP(  l.,XiHG)/4. 

1  V<*  *  *  3£MPH‘SINMbtTAH/4-  )  J+8, 

NOISE  PLUS  KA1N  ClUTTeR 

■  f“;LU,'3l:MP*H»UISE»PCN/T0KELL 
CLBAio,  =  ^t*K*SVSLF*ANAKtA.C..T0HEU*SUMSIG»RAIN/l64.0 

4  * J.261*lCoF( 2. )*RANGfc**2) 

PN01SE  =  PN01SE  ♦  CLKaIN 

TRACK! No  ACCURACY 

S1CNV  «  T EMP+bETAV 

TRKERH  *  S«JK TP C S  1GNH**2  ♦  SIGNV**2  +  ATETS  +  ATEPS  +  ATEOM 
SEARCH  ACCURACY 

fcN  =  INTF(BETAH*FR/PSCH  ♦  ,b) 

TEMP  =  SUR  TF ( tN*SJO/ ( PN01 SE+CLUTTkll 

r^iLL  *  1#^3?('51*FLAMB/<2.*PI*T0NEIL*TEMP) 

GELR  *  i.7J2o31*PUilSE*C/(2.*PI*TEMP1 

OilA  =  PNTRK*CRl*S0RTF(BETAH’#*2<-bETAV**2l/(FKS*TEMPI 

STORE  OETECTION  PARAMETERS 


CUP)  icO)  =  PNG! St 

CUP ( 99 i 

=  CLBh 

Cup ( 9  a ) 

=  BtTAH 

C OP ( 9  7 ) 

=  PULSE 

COP ( 96 ) 

=  FSCh 

CUP ( 95 ) 

=  ST  T 

CUP ( 94 ) 

sue 

CUP ( 9 j ) 

=  CLLTTk 

STLKt  OUTPUTS 

H(K,1)  = 

•  TIME 

H  t  K  » 1 1 ) 

=  FUO 

H  (  R  » 12  ) 

=  FOU-FUB 

HU,  13) 

=  oT (L ,24) *1 . E6 

H(K,  14) 

=  DtLK 

H(K,  15) 

=  UfctD 

H(K,  16) 

x  Ut  LA 

H( R,  17) 

=  SIGNV 

H(R, 16) 

*  SI  GIMH 

H(K, 19) 

=  IRKfcRk 

UU  8o  I* 

16,19 

dC  H(K,  1)  x  H(K  f I ) *1«  E3 

HU, 2c)  =  GT<c»29)*i.fc-3 
H ( K  ,  21 )  =  SGC 
Fi  U  » 22  )  *  Sub 
H ( K i 23 )  =  Sd3 
HU,  <:4 )  *  CLP 
H ( K  f  25  )  *  CLS 
HU,co)  =  CLLTTR 
OU  62  1=21, 2t 

b*  H C  K  ,  II  =  1G. *LOGlLF(H(K,  1)*1.E3) 
HU, 27)  =  lu.*LUGloF(SIG2) 

H  U  ,  28  )  =  PHA 
HU, 29)  =  PHAuLL 
HU,3o)  =  Ktf  L 


HU,j1)  *  CLBln, 

HU, 32)  *  lU.*LUGiuF<STT/(PftUISb+CLUTIR) ) 
IF  (KhT-1 )  9C,9u,S9 
9L  1PR  (  1  )  =  F LAMB 
TPk(2)  =  ANAREA 
FPROi  =  iO. *LGGluF ( AG AI NA ) 

TPR(A)  =  BETAH*5  7.2957  79  5 
TPR(5)  =  d£TAV*57. 2957795 
TPR(6)  «  FSCFl*57.2957 795 
TPk(9)  =  PI 


VPH<  10)  =  PA  VG 
TPk( 7)  =  PULSt 

IPK(il)  =  io.*LGGioF( PNOISE*!. E3 ) 

TPR ( 8 )  *  tN 

Ifl/KIlt  OUTPUT  TAPt  £>»91»(TPK(i  if  1*1,11) 

91  FURMAi  Clhl  54X  idHCOMPUTcU  CONSTANTS  ///9H  LAMBOA 
1  9h<j  Ak  =  t!3.4,oH  SOFT.  5X  bHGAIN  =  t!3.4,3H 
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GhBETaH 

9FoFSCH 

BUN 

oF  wA  I TS 
Cul\T  INOt 


x  fc  1 3  •  4 »  5H  JCG.  6X  BHBtTAV  * 
=  tl3. 4,4H  0/S  7 X  tiHTAU  « 
*  ti3.A/9HoP0  x  £13.4,7H 
5X  6HPNUISE  »  E1J.a,4H  U8M  ) 


>  £13.4, 5H 
'  E13.4.5H 
*/S«FT  4X 


»  E 1 3  •  4 , 4F) 
UB  8X 


FT./ 


DEG.  / 
SEC.  6X 
8HPAVG 


El  3.4 , 


4  „i  , 


RETURN 

END 


! 


SUEkCO T I  Nt  ECDlt 

OlFtNSiU N  Xb  AK ( O  ,  SbiML  1 2 )  , PD(2) 

OIFfcNiiUN  EPH(  322*  6)  ,CR0SEC(itiol,2),GT(b0,32),H(BOy40),  COM  l  loG  ) 
COMMON  TlMt,K3TeP,KHT,fcPH,CKOS£C,OT,H,CQM 

c 

C  GET  STUFF  FKUM  UAVL 
C 

K  *  KSTtP 
IF(NHT-i)  tiO,6u,82 
be  PuCMl  *  U« 

PULP  T  =  0. 

61  DU  b  1*1,2 

5  SohL(I)  =  CUMU  +  5i) 

CLL1TK  *  CCM(Si) 

PNUlSL  *  CCMfieO) 

UtlFC  =  U0MI99) 
ttclAh  *  CuM(Sb) 

PULSE  *  CUM  197) 

FSUH  =  CUM9t) 

C 

C  Utl  STUFF  FRCM  COMMON 

C 

Fk  *  CU.M (24) 

TACFA  *  CUM ( 3 o ) 

KASe  *  UCMtbSJ+.S 
IF(XASc)  3,3,4 
2  KASc  =  1 

4  COMINUE 

POMS  =  cUM(Se) 

L  =  K 

ALPnAb  *  GT  (  L ,  i*. ) 

C 

C  ScT  UfcTtCTlGN  PAKAMtTfckS 
C 

N  *  bcTAH*pR/FS0P  ♦  *3 
NC  *  btTAH#UfcLFC/FSuH  +  ,ti 
IF  (NO  7  » 7  » b 
7  NC  =  1 
t  CN  =  N 
cimC  =  NC 
tMe  *  cN/LNC 

FaF  *  EN*PULSfc*.69314718/TAUFA 
FAPu  =  LUGlUP(FAP) 

ZuAk  =  CLUTTK/FNC1  St 
if(FUrtS)  t>e,c,5u 

6  lFUEAR-.i)  i,o,9,S 

5  IF ( ZfcAK— 10  •  )  iu  ,  lo ,30 
C 

C  NtITPhk  uLUTTtK  NCR  NOISE  DOMINANT — fctjS.  lo-17A. 

C 

1C  CuMlNUt 

trtcZ  =  EMC*ZbAR 
VZ  *  7.*cN*ZfcAR 

PeZ  =  uOuiuF  (PFa(N»NC,Z8AR »VZ) ) 

*1  VI  *  VZ  ♦  .0*Vz*SIGNFU.,PLZ-FAPL) 

PL  1  =  L00lUP(PFA(N,NC,ZBAK,Vl)  ) 

IF ( S IGnF ( 1.,PLZ-FAPL)-SIGnF( 1,  ,FAPL~PLi) )  12,14,12 
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1*  Vi  =  Vi 
P Li  =  PL i 
GO  TC  11 

14  Vu  =  VZ  +  ( V 1- VZ  )*  (PLZ-FAPL )  / ( PLZ-PLi ) 

Du  Zi  i  =  i »  i 

X  d  A  K  l  i  )  =  SGNLd  )/P.\jUISt 
A 1  =  fc\*XuAK (I ) 

A2  =  tMC Z  *  A1 
A3  =  i.  +  A2 

b  =  c^LbOMAj)  -  LUoFiAZ)  -  (LNC- 1  •  )  *LuGF  (  A 1 )  -  VI1/A3 
Pui  1  )  =  fcXPF ( d ) 
zu  CUNT INUfc 

HlKtfc)  =  c. 

Gu  TC  99 
C 

C  CLlIItK  ULPINANT  — fcUS.  18-19 
C 

30  VZ  =  (cNC+lo.l^eKC 

PLZ  =  LOGlCF(PFAC<N,NC,VZ)  ) 

31  Vi  =  VZ  ♦  •  Gl<‘VZ<tSIGNF{i#,PLZ-FAPL) 

PL  1  =  LuGlOP  (PFAUN,NC,V1)  > 

IF  (SlG.vFil.,  PLZ-FAPL) -SI  GNP<  U,  F  APL-PL  i)  )  32,34,32 

32  VZ  =  VI 
PL  2  =  PL1 
GU  TL  31 

34  Vd  =  VZ  ♦  ( V  1-  Vi  )* (PLZ-F APL I  / 1  PLZ-PL1 ) 

DU  4L  1=1, Z 

XdAK(i)  =  'SGNLI  I  J/CLUTTR 
IF (X8Ak( I i-l.t-2)  si, 42, 42 

41  Plilll  =  c. 

GU  TC  40 

42  A 1  =  tNC*XdAR(I) 

AZ  =  (l.+Ali/Al 
A3  =  Vd/tHC 
IF(NC-l)  33,33,37 

35  PD(i)  =  fcXPF (-A3/I l.+Ai) ) 

GU  TO  HO 

37  SUP  =  0. 

JMAX  =  NC-i 

DU  3a  J= 1 , JMAX 
R  =  0  • 

DU  36  M= 1 , J 

3c  k  =  K  ♦  t VAL (A3 , M—l ) 

38  SUM  =  SUM  +  fcXPF ( FLOAT FI  NC-l-J )*LOGF ( A2I+LOGF (R)-LOGF  <  A 1 ) i 
PU(I)  =  SUM  ♦  fcXPF (-A3/I 1«  +  A1) +{ tNC— 1* ) *LUGF( A2) ) 

40  CUNT INUfc 

H(K,fc)  =  -1. 

GU  TC  99 
C 

C  NUlSfc  DOMINANT— MARCUM-SWfcRLING 
C 

SO  CUNT  INUfc 

FAN  =  LUGloFI TAUFA/I EN*PULSfc) ) 

DU  t>o  1=1,2 

Xd AK  ( 1 )  =  SGNL(I)/(CLUTTR+PNQISfc) 

CALL  MARCUM  (N,F AN, XBAR(I) ,KASE,PD(I),VBI 
60  CUNT INUfc 

H(K,t)  =  KASt 
99  CuNT INUfc 
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<-  <■>  r»  o  o 


IFIALPnAu)  199, 29*, 299 
199  POt  1  )  =  v « 

PD  (2  )  =  «,. 

299  COM  INLE 

PCLMi  =  PCUMi  ♦  ( i.-PCUMl) ♦Put  1) 
PCUil  =  PCLMT  +  (l.-PCUMT)*P0(2) 
h(R,2)  =  HU  ( 1 ) 

H(R»2)  =  PCUMl 
H  (  *  ,  4  )  *  F  0  (  2 ) 

H(K,S)  =  FCUMT 

H  (  K  f  ?  I  =  va 

H(K,b)  =  lu.*LC01CF{ SGNL (2 )/FNClSE) 
H(K,9)  *  1C.*L0' -luFIZBAK) 
h(KflO)  *  EMC 
lFIKhT—  i)  1UL  f  il'U  f  iw2 
lui  WRITE  OUTPUT  TAPE  6,lul,FAP 
101  FORMAT  ( 9HCPF A  =  Eli.4) 

1C 2  CONTINUE 
Ke  TURIN 
EnC 


4.1 

FUNCTION  PFA(N,NC,ZaARf Va) 

CCMPUTtS  FaLSE-ALARM  PRCBAbiLlTY  OF  RECEIVER  NOISE 
ANO  CLJTTEk  ONLY.  N  I b  NO.  OF  PULSES,  NC  IS  NO.  OF 
iNUtPcNUENT  CLUTTER  GROUPS,  ZBAK  IS  AVERAGE  CLUTTER  TO 
in 0 1 S  E  RATIU,  V6  IS  THE  ThKESHOLO-*  SEE  BBD-1387. 

FMC  *  FlCATF (N)/FLUATF (NC) 

FMCZ  =  FMC*ZEAR 
FHCZU  =  i.+FFCZ 
NinC  =  N-NC 
FNNC  *  FLUATF(NNC) 

SUM  *  U. 

0  »  FNNC*LOGF( FMCZU/ FMCZ )  -  SUMLUG(NNC~1 ) 

UU  99  1  =  1 ,  NC 
NU  *  1-1 
FNC  *  FLOATF(NU) 

R  *  C. 

KK  =  NC-NU 
DO  4  R=1,KK 

4  R  =  K  +  £VAL(Vd/FMCZU,K-l) 

IF (k )  99,99,5 

5  Tckrt  =  U+SUMLOGI NNC-l+NU )-SUMLOG ( NU )-FNU*LGGF (FMCZ) *LOGF ( R ) 
TERM  =  EXPFI TERM) 

I F (XMUUF ( NU, 2 ) )  1U, 20,10 

10  TERM  =  -TERM 
20  SUM  =  SUM+TERM 
99  CONTINUE 
PFA  *  SUM 
RETURN 
ENC 


4.2 

FUNCTION  PFAC ( N » NC » YC ) 

C 

C  FALSE  ALARM  PROBABILITY,  SIMPLIFIED  FOR  THE  CASE  WHERE 
C  CLLTTtR  IS  DOMINANT.  (CLO-4-117,  EQ,  18) 

C 
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4.3 

SUBKCiOUNt  MAKCUH  { N  ,F  AN  ,  SNK,KASE,  PN »  B  IAS  1 

c 

C  CUKPUTfc  MARCUM-SWtRLlNG  UfcTfcCT  I  UN  PROBABI  L I  T  I  tS 
C 

C  TEST  INPUTS 

c 

1KM  99,99,2 

2  IF(FAn)  99,99,3 

3  iFlKASfc)  99,4,4 

4  1 F (KASfc— 4 )  3,5,99 
b  1FCSNRJ  99,99,0 

C 

C  LSTiMAIc  BIAS  LCVeL 

C 

0  6  ENFh  =0* 

cN  Pf\  =  PaN 
£N  =  N 

U  YBPK  *  0. 

I F IN-12 )  7,7,B 

7  YBPK  *  fcN*(i.+2.2«ttNPK/6N**(2./i.  +  .Clb*i:NPR)  ) 

GO  ft  Li 

fc  YBPk  *  fcN*  ( 1  ,  +  i.  3<=tNPR/cN**  (  *  5*.  ol  l*fcNPk  )) 

c 

C  CUFPU1E  BIAS  LtVfcC 

C 

0  ii  fcNPh  *  lu.**ENPF 
U  UAF.Pk  =  UoAM  (  YBPK  ,  N- 1 1 

C  PYB  *  .5**(i ./tNPk) 

U  It-  IuaMPk-PYU  )  It ,12, 12 

C  1C  H  *  1 

Go  I C  14 
C  U  H  *  -.Cl 
C  14  Y  t  «  YuPK 
U  to  =  UE  V  Al  ( YU » N—  1 ) 

0  1 u  Y l  *  Yu+H 

U  ti  *  OuVAl ( Y1 , N— 1 ) 

C  SI tP  =  oAMPk  +  h*{ /2  • 

0  1F(S1GNM1.,STEP-PYB)-S1GNF(1.,H)  )  IB  ,2C  ,  lti 

C  IB  Yu  =  Yi 
u  U  =  tl 

D  GAPPK  •  SUP 

Go  1C  io 

U  2C  IF  (h)  22,24,2h 

D  22  Yb  *  Yl  -  H*(PYb-STtP)/( GAM PH- ST Er ) 

Gu  TC  it 

G  <;4  Yb  =  Yu  +  H*(  PYb-GAMPR J  /  (  STtP— QAMPK) 
it  BIAS  -  Yb 
C 

C  Stl.fcCT  P,-S  CASE 

C 

X  =  SNK 
K  *  KASt*i 

Go  1C  (ltv.,2tu,3Ut,AUu,50U)  ,  K 
C 

C  CASE  0 

c 

let  sop,  =  t. 
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P  =  tM*X 

if(Yd-p-tu)  i 5o,jc*,iu2 

IsjZ  Kb  =  -lti\!+l.  1/2.  +  bORTFU  l fcN-i. ) /2. )**2+P*YB ) 
KS  *  XMAXtF(Kb,0) 

Gb  =  1.-gAM(  Yb,KS  +  N-l ,  Ti\) 

Tb  =  tVAL(P,Kb)*Gb 
G  =  0  b 
K  =  KS 
TthM  *  Tb 
Tl  =  IN 

lie  TCPP  =  bOM+T£.KM 

If(buK-TLMP)  11^,116, iio 
112  bJP  »  TtMP 

1 f ( K )  tl6, 116,114 

114  TfcKM  *  f  tKM*f  LLATF  <K ) #  1G” TLI / ( P*G) 

G  =  G-TL 
K  =  K— 1 

TL  =  1 L*f LGATF(K+N)/Yb 
GU  1C  110 

116  TL  =  lN*Yb/PLUATHKb+N) 

K  =  KS  +  1 
G  =  Gb+T L 

TlKm  =  T  b*P*G/ (GS^fLOATFIK)  ) 

120  ItKP  *  bUM+TLKM 

lP(SUM-TfcMP)  122,1.'»J,1S0 
122  bOK  =  TtMP 

TL  =  T  L*Y0-/f  LUATf  (  K+N) 

K  =  K+l 

TfcKM  =  TfcKM#P*(G+Tlui/(G*FLOATP4KU 
G  =  G+TL 
GiJ  TL  120 

lbo  Kb  =  -1.  -  fcN/2.  +  bURTF(EN**2/4.+P*YBl 
KS  =XMAXOF (KSfO ) 

Ob  *  GAMt Yb,Kb+M-i,TN) 

IFiGSJ  174,174,155 
155  TS  *  t VAl ( P,  Kbl*GS 
G  =  Gb 
TfcKM  =  Tb 
K  =  Kb 
Tl  =  TN 

lot  TEPP  =  SUM+T fcRM 

IFlbUM-TtMPI  162 1  166  *  loo 
lo2  bUM  =  TtMP 

1 F  C  K  >  100*166, 164 

lc A  TERM  =  ieftM*FL0ATFtKi*<G+TU/<P*G4 
G  =  G+TL 

TL  =  TL*PL0ATF(K+N-i)/Y3 
K  =  K-l 
GO  TC  160 

lot  TL  =  T.Y*rii/FLGATF<K$*Hi 
K  =  Kb  +  1 
G  *  6 b~  T  L 

TfcKM  =  rb*P*G/(GS*FL£ATFUn 
170  TtMP  =  bOM  +  TERM 

IFlbLK-TfcMPI  172*174,174 
172  buP  *  TEMP 

TL  =  TL*Yb/FLUATF(K+NI 

TfcKM  *  Ti:KM*P4(G-TU/(G4FL0ATFCK+lll 

G  *  6-TL 

ri  ,  J  -*  „  \  t  $1  { \  5ft’.  5  )  ’’  t  .  Ck«-  *  •  t  >  j*r'  *  1  ■  « 


c 

c 

c 


K  =  N+  i 
Gu  f  L  i  7 1 
1/4  bUP  *  l.-bUM 

Pk  =  Sum 
GO  T  L  •ju 

bASc  i 

cut  IHiV-i)  clCfiiot22G 
tU  Phi  =  cXPF  (-YB/U.X-XI  ) 

Gu  TC  'ju 

22b  lfcPP  =  1.  +  l./<h.x*X) 

PN  =  i.  -  GAMYo,N-2#L>0m  ♦  tXPH  (fcN-i.  i*tOGF(TtMP}-Y8/U.+tN*XI  ) 
J-  *GAM(  Yb/T£ivtP»i'J~2tGUM) 

ub  Tb  4u 

C  bASt  c 

C 

3ul  Ir(N-i)  3K»31v»32u 
iiv  PiX  *  tXPF(-Yfc/(l.+X)  } 

GU  TC  So 

32u  Pim  =  i«  -  G.AM  Yb/(l.  +  XI,N-”itCuM) 

Gb  Tb  bU 
C 

C  CASE  3 

C 

4ub  IF(l\-c)  4iu»44U»43u 

4K  FiX  =  ( i.+c,*X*Yb/(X+2.)**2)*EXPP(-’2.*Y8/i2e  +  X)) 

GU  1 1  S'U 

42b  P  ix  *  (  l.+YB/ll.+X)  )*tXPF(-Y8/(  l.+XJ) 

Gb  TC  'K 

43t  C  »  2./(2.+LN*X) 
u  »  i.-u 

lF(Ye*0-lM  44u,4 5w4SU 
44u  SuP  *  U* 

ItKP  *  i. 

J  =  IS 

442  TCPP  *  SCM  +  T fckM 

IF  <  SGM-T tMP )  444  f 446*446 
444  SUP  =  ItNP 

IfckM  *  TCkM* YB^D/FCUATFl J) 

J  *  J+l 
»  ’  Gb  TG  442 

44c  Piv  *  1#  -  GAM(  Y i>  »»S-2  tUUM )  ♦  C*YB*fcVAUYB»N-2) 

1  +  u*fcVAL ( Y3»  N-l) *(  i*  +C*Y3- (fcN-2»  )  *C/D)*SUM 

GU  J C  $G 

4SC  PN  *  i.  -  GArt< YB ,N-3»0UM)  Y6*fcVAL(  Y0»N-31*C/D 

1  +  EX  PF  C-C*Yb- ( cN-2 • i ♦LUGF  {  U) )  *(  1.  ♦C<lYB‘-(  EN-2»  ) *C/0) 

2  *GAP( YB*0(N-3 »UU«) 

GU  TC  5u 

C 

b  CASt  4 

c 

Sue  SUP  =  0. 

C  *  2 • / C 2.+X) 

0  *  l.-C 
Q  »  C/U 
P  *  C«YB 

Ki  a  (3.*CN+(YB*U))  /2*~SQRTM  (  EN-l.+i  YB*D)  )**2/4.M  YB*U)*(  EN+1  .  ) ) 
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Kb  =  XMI  NcF  ( Kb , N ) 

Kb  =  AMAXoF(KS,U) 

K  a  KS 
J  =  N-Kb 
FKS  =  KS 

K  -  XMINGF (KS,N) 
if  (Y6-tN*(  l.+U))  55u,501,501 
501  Gb  =  1.  -  GAM(P,2*N-i-KS,TN) 

1 P (GS)  526,526,302 

3 02  TS  =  EXPP ( FKS*LOGF (C ) ♦( EN-FKS ) *LOGF( U) +SUMLOG( N) -SUMLUG( KS ) 
1  -SOM LOG ( J) +  LCGF ( GS) ) 

G  s  Gb 
icKM  -  TS 
TL  =  TN 

51c  TEMP  =  bUM+TERM 

IF(SOM-TEKP)  512,516, 516 
312  SUP  =  TEMP 

IF  (K )  51q,516,31A 
51A  TL  =  TL*P/FL0ATF(2*N-K) 

TERM  ^  TtRM*FLUATF  (K )  *  (G+T L )  /  <  li*FLOATF  ( N-K+l )  *G) 

0  *  G+TL 
K  *  K- 1 
GO  TO  510 

31c  IF(kS-N)  516,526,526 

316  TtRM  =  TS*w*PLUATFL<V-KS)*(G3-TN)/ ( FLUATF < KS+ 1 ) *GS) 

G  =  Gb— T  N 

TL  =  TN*FLUATF(2*N-1-KS)/P 
K  *  KS+i 

52c  TEMP  *  SUM+TERM 

i f ( som- temp )  522,520,320 

522  SUM  =  TEMP 

IF(K-N)  52a , 526 , 52o 

52 A  TERM  =  TERM*i,*FLGATF ( N-K ) * ( G-TL ) / ( FLOATF ( K+l ) *G) 

G  *  G-TL 

TL  *  TL*FLGATF(2*N-1-K)/P 
K  *  K  +  l 
Go  TO  52G 
526  PN  =  SUM 
GO  TO  50 

330  Gb  =  GAM l P , 2*N-1-KS, TN) 

IF(GS)  376,570,332 

53^  TS  =  tXPF(FKS*LUGF(C)MEN-Fl<S)*LOGF(0)«-SUMLUG(N)-SUMLUG(KS) 
i  -SUMLUG( J ) +LCGF ( GS) ) 

0  =  GS 
TERM  =  TS 
TL  =  TN 

5oG  TEMP  *  SUM+TERM 

I F ( SUM-TEMP )  562,566,560 
362  SUM  =  TLMr 

1 F ( K )  5o6,5o6,5oA 
30 A  TL  =  TL*P/FLOATF (2«N-K) 

TERM  =  TEKM*PlCATF (K)  *  (G-TL )  /  U*FlOATF  (  N-KU )  *G) 

G  =  G-TL 
K  =  K-l 
GU  TO  36U 

306  1F(KS-N)  566,376,376 

Sot  TERM  =  f  b*ii*FLUA  TF  (N-Kb)  *(  GS+TM  /( FLOATF  ( KS*1 )*G$) 
o  *  Gb+TN 

TL  =  TN*FL0ATF(2*N-1-KS) /P 
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K  =  KS  + 1 

57c  TCPP  =  SUM+fERM 

1MSUM-TLMP)  572,57o,57o 
572  SUP  =  TEMP 

IF  (K-N)  374,57t,576 

574  TERM  =  IEkH*c*FLuATMN-M*(G+TL)/(FLQATF<K+1)*G) 

G  =  G  +  TL 

TL  =  Il*FLCATFU*N-1-K)/P 

K  =  K+i 
GJ  rt  57 u 
576  PN  =  i.-SUM 
Go  Tl  sc 
C 

C  SET  PROBABILITY 

C 

SC  1  F  (PM  S  1,54,42 

91  PiM  =  c. 

GLi  TL  94 

92  lPiPiv-i.)  9*,  94,  93 
95  P.M  =  1. 

94  RETURN 
C 

L  lKRUR  MESSAGE  FGK  dAU  INPUTS 

C 

95  WRITE  OUTPUT  TaPE  o,9, N,FAN, SNR, RASE 

5  FORMAT  (lhU  /bCH  UNREASONABLE  CALL  SEQUENCE  TQ  MARCUM,  ZERO  RESULT 

1  7HS  GIVEN  //  4H  N  =  i ti,5X,5hFAN  =  ElC.8, 5X, 5HSNR  = 

2  blo«d,3X,  6HKA  SE  =  U<) 

PW  *  U. 

BIAS  =  u » 

KG  IUKiN 


n 

o 

o 

0 

0 

0 

0 

0 

0 

0 


ENC 


A  •  3  •  1 

FURL  II UN  UGaP(B,I\) 

C  INTEGRAL  =  l-( SUM,  J*C  TO  N, 

C  SUP  =  0. 

0  K  *  B 

1 F  (K-N )  lLC, 2CC, 2UC 
IOC  J  =  N+i 

C  TERM  =0EVAL( b, j  j 

C  1C  lEPP  =  SUM+TERM 
C  1 F  ( SCR- 1  EMU  15,2U,2G 

C  15  SUM  =  TEMP 
j  =  J+l 
L  EJ  =  J 

C  TERM  =  TErM*B/Fj 

GO  1L  lu 

0  Zc  Guam  =  sum 
RETUKN 
2oc  J  =  N 

C  I cRN  =  0E VaL( 6, J ) 

U  3c  ItPP  =  SUh+TERM 
C  IMSOM-TEMP)  35,40,40 

U  35  SUM  =  ItrtP 

IMJ-i)  4c , it , 3c 
0  3c  FJ  =  J 
G  TERM  *  ltkM*FJ/u 
J  =  J-l 
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-  JJ 

UF  EXPF  (  J*LUGF(  B  )-  4-LOGF  (  NFAC  )  ) 

0 

0 

0 

B 

B 

3 

3 


4.  J.  E 

FUNCTION  1/EmiVtM 
U  XPCN  =  -Y 

i F (N)  Ec i cl t  Id 
U  lo  EN  =  i\i 

C  XPCN  =  XPoN+EN*LUoF<  Y)-SUMLOG(N) 

C  Ew  UtVAL  =  cXPP  T XPON ) 

RE  I  UP  in 
eNC 


4*3*3 

PuNwTIUN  GAM(d,N,TN) 

C  SiNoLE  PRECISION  VERSION  UP  UGAM 

SUP  =  u. 

K  *  B 

I P ( K— N  )  ioC i 20C 
IOC  J  *  N+l 

TERM  =  tVAL(byJ) 

TN  *  T  ok  M*PLCATP  ( J )  /  B 
IE  TtPP  *  SUP+  TERM 

iP(SUM-TENP)  iS»<:ij»iU 
15  SUP  *  TfcMPv 
j  a  j+ 1 

PJ  =  J 

TERN  *  TtKM*6/PJ 
GO  TO  iu 
EG  GAP  *  SUM 
RETURN 
EvC  J  *  N 

TERM  *  t  V A L ( B i J ) 

IN  *  TERM 

30  TEPP  =  SUM+TEKM 

I F ( SUM-T  EMP)  35,40,40 
35  SUP  »  TEMP 

IPIJ-l)  40 , 3E, 3t 
3fc  PJ  *  J 

Term  =  term*Po/b 
J  s  J-l 
GO  TL  30 
40  GAP  *  l.-SUM 
RETURN 
ENG 


4*3.4 

FUNCTION  EVALt  Y,N) 

XPCN  =  -Y 
1 P ( N )  EG ,Eo, lo 
10  EN  =  N 

XPCN  =  XPON+EN*LOGP<  Y  Jr-SUMLOG(  N) 
EG  EVAL  *  tXPF(XPGN) 

RETURN 

ENC 


c.  c 


U  F'Ofvb  I  I LU  SOHLbb(N) 

t*  OlNtiNSlON  A(2oo) 

NM4X  =  200 

iMbOHA-UUMd  )  2b,lb»2u 
ib  OUNA  =  i. 

UUFd  =  u. 

NLAS1  *  1 
U  Ail)  *  0. 

2C  UN  =  X  ad  St  (iv  ) 

IF(NiV-i)  SLf3bf<to 

D  Jb  SUFLOG  =  b. 

Kfc  TLkN 

Hi  IF(Ms-NLaST)  bbf3u»oo 
0  5>b  SoFLbb  =  A  ( NN  ) 

KEIU'HK 

ob  K  *  MAST+i 

IF (NK-NMAX)  7b»7b*8u 
7C  00  72  i  =  K » AN 

0  72  All)  =  Atl-i)  +  LUot- (FLOATF  (  i)) 

NL  AS  I  =  NN 
GL  To  SO 

bb  l  M  I\LaST-KMAX  )  82i9ui9w 
o2  00  dA  I -K  *  UMAX 

0  bA  Ail)  =  A(i-i)  «■  LuGF (FLOATF U )) 
NLA  Si  =  UMAX 
C  9b  b  =  A(NMAX) 

K  *  MAX  +  x 

0(j  92  I*K»AN 

92  G  =  b  +  LLGF (FLuATFi I ) ) 

SOFLLG  =  b 
Kt  T0ki\ 
fcNL' 


0 

0 

0 

0 

0 

0 

0 

D 

0 

0 

0 

0 
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